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^ Abstract 



We give a pedagogical introduction to the basic ideas and concepts of the Asymptotic 
Safety program in Quantum Einstein Gravity. Using the continuum approach based upon 
the effective average action, we summarize the state of the art of the field with a particular 
focus on the evidence supporting the existence of the non-trivial renormalization group 
fixed point at the heart of the construction. As an application, the multifractal structure of 
the emerging space-times is discussed in detail. In particular, we compare the continuum 
prediction for their spectral dimension with Monte Carlo data from the Causal Dynamical 
Triangulation approach. 



*To appear in the New Journal of Physics special issue on Quantum Einstein Gravity. 



1 Introduction 



Finding a consistent and fundamental quantum theory for gravity is still one of the most 
challenging open problems in theoretical high energy physics to date [1] . As is well known, 
the perturbative quantization of the classical description for gravity, General Relativity, 
results in a non-rcnormalizable quantum theory [2-4] . One possible lesson drawn from this 
result may assert that gravity constitutes an effective field theory valid at low energies, 
whose UV completion requires the introduction of new degrees of freedom and symmetries. 
This is the path followed, e.g., by string theory. In a less radical approach, one retains 
the fields and symmetries known from General Relativity and conjectures that gravity 
constitutes a fundamental theory at the non-perturbative level. One proposal along this 
line is the Asymptotic Safety scenario [5,6] which, motivated by gravity in 2 + e dimensions 
[7,8], was initially put forward by Weinberg [9,10]. The key ingredient in this construction 
is a non-Gaussian fixed point (NGFP) of the gravitational renormalization group (RG) 
flow, which controls the behavior of the theory at high energies and renders physical 
quantities safe from unphysical divergences. 

The primary tool for investigating this scenario is the functional renormalization group 
equation (FRGE) for gravity [11], which constitutes the spring-board for the detailed anal- 
ysis of the gravitational RG flow at the non-perturbative level [11-36].^ The FRGE defines 
a Wilsonian RG fiow on a theory space which consists of all diffeomorphism invariant func- 
tionals of the metric g/^i^ and yielded substantial evidence for the existence and predictivity 
of the NGFP underlying the Asymptotic Safety conjecture. The theory emerging from this 
construction. Quantum Einstein Gravity (henceforth denoted "QEG"), defines a consis- 
tent and predictive quantum theory for gravity within the framework of quantum field 
theory. We stress that QEG is not a quantization of classical General Relativity: its bare 
action corresponds to a non-trivial fixed point of the RG flow and is a prediction therefore. 

The approach of [11] employs the effective average action [38—41] which has crucial 

^Independent support for the Asymptotic Safety conjecture comes from a 2-dimensional symmetry 
reduction of the gravitational path-integral [37]. 
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advantages as compared to other continuum implementations of the Wilsonian RG flow 
[42]. In particular, the RG scale dependence of Ffe is governed by the FRGE [38] 
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(1.1) 



Here ^ is the collection of all dynamical fields considered and $ denotes their back- 
ground counterparts. Moreover 7?.j. is a matrix- valued infrared cutoff, which provides 
a A;-dependent mass-term for fluctuations with momenta <C A;^, while vanishing for 
S> k"^. Solutions of the flow equation give rise to families of effective field theories 
{rfe [^//i^] , < k < go} labeled by the coarse graining scale k. The latter property opens 
the door to a rather direct extraction of physical information from the RG flow, at least 
in single-scale cases: If the physical process under consideration involves a single typical 
momentum scale po only, it can be described by a tree-level evaluation of T^lg^i,], with 



A striking consequence of the scale-dependence in is the observation that the ef- 
fective QEG space-time should have certain features in common with a fractal [13,15]. 
The property underlying this assertion is that the effective field equations derived from 
the gravitational average action equip every given smooth space-time manifold with, in 
principle, infinitely many different (pseudo) Riemannian structures, one for each coarse 
graining scale [43,44]. Thus, very much like in the famous example of the coast line of Eng- 
land [45] , the proper length on a QEG space-time depends on the "length of the yardstick" 
used to measure it. Earlier on similar fractal properties had already been found in other 
quantum gravity theories, in particular near dimension 2 [46], in a non-asymptotically safe 
model [47], and by analyzing the conformal anomaly [48]. 

In ref. [13] the consequences of this scale-dependence for the 4-dimensional graviton 

propagator has been studied in the regime of asymptotically large momenta and it has 

been found that near the Planck scale a kind of dynamical dimensional reduction oc- 

^The precision which can be achieved by this effective field theory description depends on the size of 
the fluctuations relative to mean values. If they turn out large, or if more than one scale is involved, it 
might be necessary to go beyond the tree-level analysis. 
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curs. As a consequence of the NGFP controlling the UV behavior of the theory, the 4- 
dimensional graviton propagator essentially behaves 2-dimensional on microscopic scales. 
Subsequently, the "finger prints" of the NGFP on the fabric of the effective QEG space- 
times have been discussed in [15], where it was shown that Asymptotic Safety induces a 
characteristic self-similarity of space-time on length-scales below the Planck length ipi^. 
In particular the spectral dimension of the effective QEG space-times in the asymptotic 
scaling regime implied by the NGFP becomes dg = d/2 with d the classical space-time di- 
mension [13,15,49]. Based on this observation it was argued in a cosmological context that 
the geometry fluctuations originating from the scaling regime can give rise to a scale free 
spectrum of primordial density perturbations responsible for structure formation [50,51]. 

Along a different line of investigations, the Causal Dynamical Triangulation (CDT) 
approach has been developed and first Monte Carlo simulations were performed [52-57], 
see [58] for a recent review. In this framework one attempts to compute quantum gravity 
partition functions by numerically constructing the continuum limit of an appropriate 
statistical mechanics system. From the perspective of the latter, this limit amounts to 
a second order phase transition. If CDT and its counterpart QEG, formulated in the 
continuum by means of the average action, belong to the same universality class^ one may 
expect that the phase transition of the former is described by the non-trivial fixed point 
underlying the Asymptotic Safety of the latter. 

Remarkably, ref. [54] reported results which indicated that the 4-dimensional CDT 
space-times, too, undergo a dimensional reduction from four to two dimensions as one 

"zooms" in on short distances. In particular it had been demonstrated that the spectral 
dimension dg measured in the CDT simulations has the very same limiting behaviors, 
4 — )• 2, as in QEG. Therefore it was plausible to assume that both approaches indeed 
"see" the same continuum physics. 

This interpretation became problematic, however, when it turned out that the Monte 
Carlo data corresponds to a regime where the cutoff length inherent in the triangulations is 

^For the time being this is merely a conjecture, of course, albeit a very natural one. 
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still significantly larger than the Planck length. The situation became even more puzzling 
when ref. [56] carried out CDT simulations for d = 3 macroscopic dimensions, which favor a 
value near = 2 on the shortest length-scale probed. Furthermore, the authors of ref. [59] 

reported simulations within the euclidean dynamical triangulation (EDT) approach in 
d = 4, where the spectral dimension dropped from 4 to about 1.5. Obviously, both of 
these observations conflict the QEG expectations if one interprets the latter dimension as 
the value in the continuum limit. 

In order to resolve this puzzle, [60] computed several types of scale dependent effective 
dimensions, specifically the spectral dimension dg and the walk dimension dy^ for the 
effective QEG space-times. Surprisingly, the analysis revealed a further regime which 
exhibits the phenomenon of dynamical dimensional reduction on length scales slightly 
larger than £pi^. There the spectral dimension is even smaller than near the fixed point, 
namely ds = 4/3 in the case of 4 dimensions classically. Moreover it was found that the 
(3-dimensional) results reported in [56] are in good agreement with QEG. This analysis 
also confirmed the supposition [56] that the shortest possible length scale achieved in the 
simulations is not yet close to the Planck length. Rather the Monte Carlo data probes the 
transition between the classical and the newly discovered "semi-classical" regime. 

It is intriguing that Loop Quantum Gravity and spin foam models also show indi- 
cations for a similar dimensional reduction [61-63], with some hints for an intermittent 
regime where the spectral dimension is smaller than in the deep ultraviolet. In ref. [64] an 
argument based upon the strong coupling limit of the Wheeler-DeWitt equation was put 
forward as a possible explanation of this dimensional reduction. Within non-commutative 
geometry Connes et al. [65] interpreted the dynamical dimensional reduction to dg = 2, 
which was observed in QEG, in the context of the derivation of the Standard Model from 
a spectral triple. In fact, from the data encoded in a spectral triple, its Dirac operator in 
particular, one can compute a type of spectral dimension of the resulting non-commutative 
space which is closely related to the one considered here. Also for standard fractals such as 
Cantor sets, it has been possible to find spectral triples representing them and to compute 
the corresponding dimensions [66] . 
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Furthermore, a number of model systems (quantum sphere, K-Minkowski space, etc.) 
give rise to a similar reduction as fully fledged quantum gravity [67] . Among other devel- 
opments, these findings also motivated the investigation of physics on prescribed fractal 
space-times. In refs. [68, 69] a fractional differential calculus [70] was employed in order 
to incorporate fractal features, and in [71] recent exact results on spectral zeta-functions 
on certain fractals [72] were used to study the thermodynamics of photons on fractals. 
In ref. [73] matter quantum field theories were constructed and renormalized on a fractal 
background. The almost universal appearance of fractional properties of space-time and 
its accessibility in various, a priori different, approaches to quantum gravity make the 
generalized notions of dimensionality which we are going to review here a valuable tool in 
comparing the physics content of these different formulations. 

The present article is intended to provide the necessary background for understand- 
ing the frontier developments in Asymptotic Safety. In the next section we follow [5] 
and introduce the general concepts related to the Wilsonian picture of renormalization: 
theory space, renormalization group flows, non-perturbative renormalizability, and the 
non-perturbative approximation scheme of truncating theory space. The construction of 
the effective average action and its FRGE for gravity [11] is reviewed in section 3. In sec- 
tion 4 we illustrate the most commonly used approximation scheme of "truncating theory 
space" by means of a simple example, the so-called Einstein-Hilbert truncation. Section 5 
contains a state of the art summary of the results obtained using truncated flow equations, 
with an emphasis on the question as to whether there exists a non-trivial fixed point for 
the average action's RG flow. If so, QEG could be established as a fundamental theory 
of quantum gravity which is non-perturbatively renormalizable and "asymptotically safe" 
from unphysical divergences. The remainder of the review is dedicated to the discussion of 
the multifractal space-times emerging from QEG: section 6 reviews the structures under- 
lying these fractal properties, the scale-dependent space-time metrics. Subsequently, the 
spectral, walk, and Hausdorff dimensions introduced in appendix A are computed in the 
framework of QEG in sections 7 and 8. In particular subsection 8.3 compares the QEG 
results for the spectral dimension with the one found in the CDT program. We close with 
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a short summary and some concluding remarks in section 9. 



2 RG flows, theory space, and Asymptotic Safety 

The key idea of the Wilsonian renormahzation group is the description of a physical system 
in terms of a one-parameter family of effective actions, each valid at a certain energy scale 
k. The RG flow of the theory thereby connects the effective descriptions at different 
scales. The arena in which this dynamics takes place is the "theory space" sketched in 
fig. 1. In order to describe it, we first specify our "theory" by fixing its field content 
^(x) and, possibly, imposing certain symmetry requirements (a Z2-symmetry for a single 
scalar, or diffeomorphism invariance if $ denotes the space-time metric, for instance). The 
theory space corresponding to this theory consists of all (action) functionals A : $ 
depending on this set of fields and compatible with the symmetry requirements. Thus the 
theory space {^[ • ]} is fixed once the field content and the symmetries are given. 

Let us assume we can find a set of "basis functionals" {-Pa[']} so that every point of 
theory space has an expansion of the form 



The basis {-Pa[']} will include both local field monomials and non-local invariants and 
we may use the "generalized couplings" {ua,a = 1,2, •••} as local coordinates. More 
precisely, the theory space is coordinatized by the subset of "essential couplings", i.e., 
those coordinates which cannot be absorbed by a field reparameterization. 

Geometrically speaking the FRGE for the effective average action, cq. (1.1), defines 
a vector field (3 on theory space. The integral curves along this vector field are the "RG 
trajectories" k parameterized by the scale k. They start, for A; ^ oo, at the 

microscopic action S and terminate at the ordinary effective action at A; = 0. The natural 
orientation of the trajectories is from higher to lower scales k, the direction of increasing 



oo 




(2.1) 



7 




Figure 1: The points of theory space are action functional A[-]. The RG equation defines a 
vector field (3 on this space; its integral curves are the RG trajectories fc i— > r^. They emanate 
from the fixed point action F, and end at the standard effective action F. 

"coarse graining". Expanding as in (2.1), 

oo 

rk[<i>,^ = Y.Ua{k)Pa[<^,^, (2.2) 
a=l 

the trajectory is described by infinitely many "running couplings" Ua{k). Inserting (2.2) 
into the FRGE we obtain a system of infinitely many coupled differential equations for 
the Uq's: 

kdkUaik) = Paiui,U2,- ■ ■ ;k) , a = l,2,--- . (2.3) 

Here the "beta functions" (3^ arise by expanding the trace on the RHS of the FRGE in 
terms of • ]}, i.e., ^Tr[---] = Yl'^=i (^a{'^ij^2, ■ ' ' ',k)Pa[^,^]- The expansion coeffi- 
cients /3„ have the interpretation of beta functions similar to those of perturbation theory, 
but not restricted to relevant couplings. In standard field theory jargon one would refer to 
Uaik = oo) as the "bare" parameters and to Ua{k = 0) as the "renormalized" or "dressed" 
parameters. 
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The notation with the bar on and is to indicate that we are still dealing with 
dimensionful couplings. Usually the flow equation is reexpressed in terms of the dimen- 
sionless couplings 

Ua = k-'^'^ Ua , (2.4) 

where da is the canonical mass dimension of Uq. Correspondingly the essential n^'s are 
used as coordinates of theory space. The resulting RG equations 

kdkUaik) = Paiui,U2,- ■ ■) (2-5) 

are a coupled system of autonomous differential equations. The /3q's have no explicit 
/c-depcndence and define a "time independent" vector field on theory space. The RG 
trajectories arise as solutions or, equivalently, integral curves of (2.5). 

Based on these structures, the concept of renormalization can be understood as fol- 
lows. The boundary of theory space depicted in fig. 1 is meant to separate points with 
coordinates {ua,a = 1, 2, • • • } with all the essential couplings Ua well defined, from points 
with undefined, divergent couplings. The basic task of renormalization theory consists 
in constructing an "infinitely long" RG trajectory which lies entirely within this theory 
space, i.e., a trajectory which neither leaves theory space (that is, develops divergences) 
in the UV limit A; — >■ oo nor in the IR limit A; — >■ 0. Every such trajectory defines one 
possible quantum theory. 

The key idea of Asymptotic Safety is to perform the UV limit — > oo at a fixed 
point {n*,a = 1,2, • • • } = u* of the RG flow. The fixed point is a zero of the vector 
field f3 = (/3q), i.e., f3a{u*) = for all a = 1,2,---. The RG trajectories have a low 
"velocity" near a flxed point because the /^q's are small there and directly at the fixed 
point the running stops completely. As a result, one can "use up" an infinite amount of 
RG time near/at the fixed point if one bases the quantum theory on a trajectory which 
runs into a fixed point for A; ^ oo. The fact, that, in the UV limit the trajectory ends 
at a fixed point, an "inner point" of theory space giving rise to a well behaved action 
functional, ensures that the trajectory does not escape from theory space, i.e., does not 
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Figure 2: Schematic picture of the UV critical hypersurface <Suv of the NGFP. It is spanned by 
RG trajectories emanating from the NGFP as the RG scale k is lowered. Trajectories not in the 
surface are attracted towards 5uv as k decreases. (The arrows point in the direction of decreasing 
k, from the "UV" to the "IR".) 

develop pathological properties such as divergent coTiplings. For A; — )• oo the resulting 
quantum theory is "asymptotically safe" from unphysical divergences. 

In a sense, standard perturbation theory takes the k — > oo-Iimit at the Gaussian fixed 

point (GFP), i.e., a fixed point where n* = 0, Va = 1, 2, This construction is the one 

underlying asymptotic freedom. More general, however, one can also use a non-Gaussian 
fixed point (NGFP) for letting k — t- oo, where, by definition, not all of the coordinates 
n* vanish. In the context of gravity, Weinberg [9] proposed that the UV completion of 
the theory is precisely given by such a NGFP, which therefore constitutes the essential 
ingredient in the Asymptotic Safety program. 

Given a NGFP, an important concept is its UV critical hypersurface S\jy, or synony- 
mously, its unstable manifold. By definition, it consists of all points of theory space which 
are pulled into the NGFP by the inverse RG flow, i.e., for increasing k. Its dimensionality 
dim (5uv) = ^uv is given by the number of attractive (for increasing cutoff k) directions 
in the space of couplings. 

For the RG equations (2.5), the linearized flow near the fixed point is governed by the 
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Jacobi matrix B = (-B^-y), B^-y = d^Paiu*): 

k dk Uaik) = ^ (u-yik) - u*) . (2.6) 

7 

The general solution to this equation reads 

Ua{k) = ul + ^CiV^ (^^y (2.7) 

where the V^^s are the right-eigenvectors of B with eigenvalues —6j, i.e., = 
—OjV^. Since B is not symmetric in general the 6/'s are not guaranteed to be real. We 
assume that the eigenvectors form a complete system though. Furthermore, ko is a fixed 
reference scale, and the C/'s are constants of integration. The quantities Oj are referred to 
as critical exponents since when the renormalization group is applied to critical phenomena 
(second order phase transitions) the traditionally defined critical exponents are related to 
the 6'/'s in a simple way [40]. 

If Ua{k) is to describe a trajectory in Suy, Ua{k) must approach -u* in the limit A; — >■ oo 
and therefore we must set Cj = for all I with Re ^7 < 0. Hence the dimensionality 
Auv equals the number of B-eigenvalues with a negative real part, i.e., the number of 
^/'s with KeOj > 0. The corresponding eigenvectors span the tangent space to 5uv 
at the NGFP. If we lower the cutoff for a generic trajectory with all Cj nonzero, only 
Auv "relevant" parameters corresponding to the eigendirections tangent to Sjjy grow 
(Re ^7 > 0), while the remaining "irrelevant" couplings pertaining to the eigendirections 
normal to iSuv decrease (Re 6j < 0). Thus near the NGFP a generic trajectory is attracted 
towards Suy, see fig. 2. 

Coming back to the Asymptotic Safety construction, let us now use this fixed point in 
order to take the limit k — t- 00. The trajectories which define an infinite cutoff limit are 
special in the sense that all irrelevant couplings are set to zero: C/ = if HeOj < 0. These 
conditions place the trajectory exactly on Syy. There is a Auv-parameter family of such 
trajectories, and the experiment must decide which one is realized in Nature. Therefore 
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the predictive power of the theory increases with decreasing dimensionaUty of Sjjy, i.e., 
number of UV attractive eigendirections of the NGFP. If Auv < oo, the quantum field 
theory thus constructed is comparable to and as predictive as a perturbatively renormal- 

izable model with Auv "renormalizable couplings", i.e., couplings relevant at the GFP, 
see [5] for a more detailed discussion. 

Up to this point our discussion did not involve any approximation. In practice, how- 
ever, it is usually impossible to find exact solutions to the flow equation. As a way out, 
one could evaluate the trace on the RHS of the FRGE by expanding it with respect to 
some small coupling constant, for instance, thus recovering the familiar perturbative beta 
functions. A more interesting option which gives rise to non-perturbative approximate 
solutions is to truncate the theory space {^[-j}. The basic idea is to project the RG flow 
onto a finite dimensional subspace of theory space. The subspace should be chosen in such 
a way that the projected flow encapsulates the essential physical features of the exact flow 
on the full space. 

Concretely the projection onto a truncation subspace is performed as follows. One 
makes an ansatz of the form r,t[^)^] = 'l2i^i ^ii.^) ^ii^ ^ ^] ^ where the fc-independent 
functional {Pi['])^ = Ir'' j-^} form a 'basis' on the subspace selected. For a scalar 
field (f), say, examples include pure potential terms J d^x<p'^{x), J ln^^(a:), ••• , 

a standard kinetic term jd'^x{d4>Y, higher order derivative terms Jd'^x ^ (5^)" ^, 
and non-local terms like / d'^x (l)\n{—d'^)<^, • • • . Even if 5 = Too is simple, a standard 0^ 
action, say, the evolution from k = oo downwards will generate such terms, a priori only 
constrained by symmetry requirements. The difficult task in practical RG applications 
consists in selecting a set of Pj's which, on the one hand, is generic enough to allow for a 
sufficiently precise description of the physics one is interested in, and which, on the other 
hand, is small enough to be computationally manageable. 

The projected RG flow is described by a set of ordinary (if N < oo) differential 
equations for the couplings Ui{k). They arise as follows. Let us assume we expand the 
^-dependence of ^Tr[- • • ] (with the ansatz for rfc[$, $] inserted) in a basis {-Pa[ ■ ]} of the 
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full theory space which contains the Pj's spanning the truncated space as a subset: 

^ oo A'^ 

-TV[---] = ^^„(«i,--- ,UN;k)Pa[^,^ = J2Mui,--- ,UN;k)Pi[<P,^+icest. (2.8) 

a=l 1=1 

Here the "rest" contains all terms outside the truncated theory space; the approximation 
consists in neglecting precisely those terms. Thus, equating (2.8) to the LHS of the flow 
equation, dtV^ = X^^^i dtUi{k)Pi, the linear independence of the P^'s implies the coupled 
system of ordinary differential equations 

dtUi{k) = Pi{ui,--- ,UN;k) , i = l,---,N. (2.9) 

Solving (2.9) one obtains an approximation to the exact RG trajectory projected onto the 
chosen subspace. Note that this approximate trajectory does, in general, not coincide with 
the projection of the exact trajectory, but if the subspace is well chosen, it will not be very 
different from it. In fact, the most non-trivial problem in using truncated flow equations 
is to find and justify a truncation subspace which should be as low dimensional as possible 
to make the calculations feasible, but at the same time large enough to describe at least 
qualitatively the essential physics. We shall return to the issue of testing the quality of a 
given truncation later on. 

3 The effective average action for gravity 

In the case of QEG, ideally, we would like theory space to consist of functionals ^[g'lijy] 
depending on a symmetric tensor field, the metric, in a diffeomorphism invariant way. 
Given a theory space, the form of the FRGE and, as a result, the vector field P are 
completely fixed. However, in the case of gravity it is much harder to make this idea 
work in a concrete way as compared to a simple matter field theory on a non-dynamical 
space-time, for instance. The reasons are of both conceptual and technical nature: 

(1) The theory of quantum gravity we are aiming at should be formulated in a background 
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independent way. It should explain rather than presuppose the existence and the proper- 
ties of space-time. Hence no special space-time manifold, in particular no special causal, 
let alone Riemannian, structure should play a distinguished role at the fundamental level 
of the theory. In particular it should also cover "exotic" phases in which the metric is de- 
generate or has no expectation value at all. Since almost our entire repertoire of quantum 
field theory methods applies only to the case of an externally given space-time manifold, 
usually Minkowski space, this is a severe conceptual problem. In fact, it is the central 
challenge for basically all approaches to quantum gravity, in one guise or another [1]. 

(2) The second difficulty, while less deep, is important for the practical applicability of 
the RG methods to be developed. It occurs already in the standard functional integral 
quantization of gauge or gravity theories, and is familiar from Yang-Mills theories. If 
one gauge-fixes the functional integral with an ordinary (covariant) gauge fixing condition 
like df^A*^ = 0, couples the (non-abelian) gauge field to a source, and constructs the 
ordinary effective action, the resulting functional r[>l^] is not invariant under the gauge 
transformations of A^, A'^ ^ A^j^ + D'^{A)uJ'. Only at the level of physical quantities 
constructed from r[74^], S-matrix elements for instance, gauge invariance is recovered. 

(3) A more profound problem is related to the fact that in a gauge theory a "coarse 
graining" based on a naive Fourier decomposition of A'^{x) is not gauge covariant and 
hence not physical. In fact, if one were to gauge transform a slowly varying A'^{x) using a 
parameter function uj"'{x) with a fast x- variation, a gauge field with a fast a;- variation would 
arise which, however, still describes the same physics. In a non-gauge theory the coarse 
graining is performed by expanding the field in terms of eigenfunctions of the (positive) 
operator —d"^ and declaring its cigcnmodcs 'long' or 'short' wavelength depending on 
whether the corresponding eigenvahie is smaller or larger than a given fe^. In a gauge 
theory the best one can do in installing this procedure is to expand with respect to the 
covariant Laplacian or a similar operator, and then organize the modes according to the 
size of their eigenvalues. While gauge covariant, this approach sacrifices to some extent 
the intuition of a Fourier coarse graining in terms of slow and fast modes. Analogous 
remarks apply to theories of gravity covariant under general coordinate transformations. 
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The key idea which led to a solution of all three problems was the use of the background 
field method [74] . At first sight it seems to be a contradiction in terms to use background 
fields in order to achieve "background independence". However, actually it is not, since 
the background metric introduced, g^^, is kept completely arbitrary, and no physics ever 
may depend on it. In fact, it becomes a second argument of [51^;^ , ^^j^] and may be 
freely chosen from the same function space the dynamical metric "lives" in. The use of 
a background field also opens the door for complying with the requirement of a gauge 
invariant effective action. It is well known [75,76] that there exist special gauge choices, 
the so-called background gauge fixing conditions, that make the (ordinary) effective action 
a gauge or diffeomorphism invariant functional of its arguments (including the background 
fields!). As it turned out [11,39] this technique also lends itself for implementing a covariant 
IR cutoff, and it is at the core of the effective average action for Yang-Mills theories [39,41] 
and for gravity [11]. In the following we briefly review the effective average action for 
gravity which has been introduced in ref. [11]. 

The ultimate goal is to give meaning to an integral over 'all' metrics 7^;^ of the form 
jT^lfiu exp{— 5[7^jy] + source terms} whose bare action «S'[7^i/] is invariant under general 
coordinate transformations. 



where is the Lie derivative with respect to the vector field v^d^. To start with we con- 
sider 7^1, to be a Riemannian metric and assume that Sl^jf^^] is positive definite. Heading 
towards the background field formalism, the first step consists in decomposing the variable 
of integration according to j^^, = g^^, + h^u, where ^^i^ is a fixed background metric. Note 
that we are not implying a perturbative expansion here, h^i^ is not supposed to be small 
in any sense. After the background split the measure V^^i, becomes Vh^^ and the gauge 
transformations which we have to gauge-fix read 



(3.1) 



5h 




0. 



(3.2) 
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Picking an a priori arbitrary gauge fixing condition Fn{h\g) = the Faddeev-Popov trick 
can be applied straightforwardly [75]. Upon including an IR cutoff Afc5'[/i, C, (7; 5] (cfg. 
eq. (3.13) below) we are led to the following fc-dependent generating functional Wk for the 
connected Green functions: 

exp C7^ a^- g^,]] = j Vh^.VC'^VC^ exp ^-S[g + h]- Sgf[h; g] 

- 5gh[/i, C, C; g] - AkS[h, C, C; g] - ^source} • (3.3) 

Here S'gf denotes the gauge fixing term 

Sdh;9] = ^J <fx^r''F^F, , (3.4) 
and ^gh is the action for the corresponding Faddeev-Popov ghosts and C^: 

/- f)T^ 
d^x^g TTT-^ (gap + Kp) . (3.5) 

The Faddeev-Popov action S'gh is obtained along the same lines as in Yang-Mills theory: 
one applies a gauge transformation (3.2) to and replaces the parameters by the ghost 
field C^. The integral over and exponentiates the Faddeev-Popov determinant 
dei[5Fn/6v'^]. In (3.3) we coupled /i^^, and 0^ to sources t^^, and a^, respectively: 
'S'source = — J df^x "v/^ji'^'^/i^i/ + d^f^C^ + cr'^C^j . The k- and source-dependent expectation 
values of h^i,, and (7^ are then given by 

As usual we assume that one can invert the relations (3.6) and solve for the sources 
{f^" , a^^ , a^) as functionals of (Ji^^ , ^'^ , C/i) and, parametrically, of g^i,. The Legendre 
transform of Wk reads 

f fe[/i, e"; 9]= J d'^x [t^'Ku + ^,.e + <^%} - Wk[t, a, a; g] . (3.7) 
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This functional inherits a parametric 5^,/-dependence from Wk- 

As mentioned earlier for a generic gauge fixing condition the Legendre transform (3.7) 
is not a diffeomorphism invariant functional of its arguments since the gauge breaking 
under the functional integral is communicated to via the sources. While Fk does 
indeed describe the correct 'on-shell' physics satisfying all constraints coming from BRST 
invariance, it is not invariant off-shell [75,76]. The situation is different for the class of 
gauge fixing conditions of the background type. While - as any gauge fixing condition 
must - they break the invariance under (3.2) they are chosen to be invariant under the 
so-called background gauge transformations 

^h/j,^ = Cyhfj^i, , 6g^y = CyQ^j^i, . (3.8) 

The complete metric = g^i, + h^i, transforms as 6j^i, = C^j^i, both under (3.8) and 
under (3.2). The crucial difference is that the ('quantum') gauge transformations (3.2) 
keep §1^1^ unchanged so that the entire change of 7^,^ is ascribed to h^i,. This is the point 
of view one adopts in a standard perturbative calculation around flat space where one 
fixes gfj,i, = 77^,/ and allows for no variation of the background. In the present construction, 
instead, we leave g^i, unspecified but insist on covariance under (3.8). This will lead to a 
completely background covariant formulation. 

Clearly there exist many possible gauge fixing terms Sgi[h;g] of the form (3.4) which 
break (3.2) and are invariant under (3.8). A convenient choice which has been employed 
in practical calculations is the one-parameter family of gauge conditions 

F^[h;g] = V2K{D''h^,-wD^h\), (3.9) 

parameterized by w. The covariant derivative involves the Christoffel symbols Tfn, of 
the background metric. Note that (3.9) is linear in the quantum field hap- For w = ^, 
(3.9) reduces to the background version of the harmonic coordinate condition [75]: on 
a flat background with g^i, = -q^i, the condition = becomes the familiar harmonic 
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coordinate condition, d^h^u = ^d,^h/^. In eqs. (3.9) and (3.5) k is an arbitrary constant 
with the dimension of a mass. We shall set k = {32TrG)~^^^ with G a constant reference 
value of Newton's constant. The ghost action for the gauge condition (3.9) reads 

Sgh[h,C,C;g] = j (X^x ^C^M[g,gT^C^ (3.10) 

with the Faddeev-Popov operator 

M [g, gf, = DP g^. Dp + D" g^.D^^ -2wD'' gf"' gp,D„ . (3.11) 

It will prove crucial that for every background- type choice of -F^, is invariant under 
(3.8) together with 

dC^^ = L^C^ , 5C^ = C^C^ . (3.12) 

The essential piece in eq. (3.3) is the IR cutoff for the gravitational field h^n and for 
the ghosts. It is taken to be of the form 

^kS = ^ j d''xVdK.nf^"[9r"'hp. + V2J d'^xV^c^nfiaic'^ . (3.13) 

The cutoff operators Uf^^ and Uf serve the purpose of discriminating between high- 
momentum and low-momentum modes. Eigenmodes of —D^ with eigenvalues S> are 
integrated out without any suppression whereas modes with small eigenvalues ^ are 
suppressed. The operators TZf^^ and TZf^ have the structure TZkld] = Z^k'^ R^^\— D'^ / k'^) , 
where the dimensionless function R^'^^ interpolates between R^^\0) = 1 and R^^\oo) = 0. 
A convenient choice is, e.g., the exponential cutoff R^^\w) = w[exp{w) — or the 
optimized cutoff R^'^\w) = (1 — w)9{l — w), where w = p^/k"^. The factors are 
different for the graviton and the ghost cutoff. They are determined by the condition that 
all fluctuation spectra are cut off at precisely the same k'^, such that TZk combines with 
r^^^ to the inverse propagator T^^^ + TZk = Zk{p^ + k'^) + ■ ■ ■ , as it is necessary if the IR 
cutoff is to give rise to a (mass)^ of size k'^ rather than [Z^j'^k"^. Following this condition 
the ghost Zk = zf^ is a pure number, whereas for the metric fluctuation Z^ = Zf^^ is a 
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tensor, constructed only from the background metric g^n. 

A feature of A^S" which is essential from a practical point of view is that the modes 
of h^y and the ghosts are organized according to their eigenvalues with respect to the 
background Laplace operator D"^ = g'^'^D^j^Di, rather than = gi^^Dfj,Di,, which would 
pertain to the full quantum metric g^i, + h^u. Using D^, the functional AjtS is quadratic 
in the quantum field hfj^i^ while it becomes extremely complicated if D'^ is used instead. 
The virtue of a quadratic AjfeS" is that it gives rise to a flow equation which contains only 
second functional derivatives of Ffc but no higher ones. The flow equations resulting from 
the cutoff operator are prohibitively complicated and can hardly be used for practical 
computations. A second property of A^S* which is crucial for our purposes is that it is 
invariant under the background gauge transformations (3.8) with (3.13). 

Having specified all the ingredients which enter the functional integral (3.3) for the 
generating functional we can write down the final definition of the effective average 
action r^. It is obtained from the Legendre transform Tfc by subtracting the cutoff action 
AfeS" with the classical fields inserted: 

Vk[ki,i\9]=^k[Ki,i\g]-^kS[h,i,^-g]- (3.14) 
It is convenient to define the expectation value of the quantum metric 7^,/, 

9,xy{x) = g^,^{x) + h^,^{x), (3.15) 
and consider Tjt as a functional of g^i, rather than h^n,'- 

rit[5/ui/, giiv, C Cm] = '^kigfj.u - gij.v, C 9tJ-v] ■ (3-i6) 

So, what did we gain going through this seemingly complicated background field con- 
struction, eventually ending up with an action functional which depends on two metrics 
even? The main advantage of this setting is that the corresponding functional , and as 
a result T^, are invariant under general coordinate transformations where all its arguments 
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transform as tensors of the corresponding rank: 



Fkl^ + Cy^] = , * = {9n^,9^iu,e,^^^} ■ (3.17) 

Note that in (3.17), contrary to the "quantum gauge transformation" (3.2), also the back- 
ground metric transforms as an ordinary tensor field: Sg^i^ = Cvg^v Eq. (3.17) is a 
consequence of 

Wk [J + C^J] = Wk[J], J = {^'^^ C7^ a^; g^,] . (3.18) 

This invariance property follows from (3.3) if one performs a compensating transformation 
(3.8), (3.13) on the integration variables /i^^, and and uses the invariance of S[g + 
h], Sg{, Sgh and A^S. At this point we assume that the functional measure in (3.3) is 
diffeomorphism invariant. 

Since the T^fc's vanish for k = 0, the limit A; — )■ of Fk[gij.iy,g^u,^^,^fi] brings us back 
to the standard effective action functional which still depends on two metrics, though. 
The "ordinary" effective action r[g'^jy] with one metric argument is obtained from this 
functional by setting g^i^ = ^^j^, or equivalently h^iy = [75,76]: 

T[g] = lim Tkig, g = g,^ = 0,^ = 0] = lim Tkih = 0,^ = 0,^ = 0;g = g] . (3.19) 

This equation brings about the "magic property" of the background field formalism: a 
priori the IPI n-point functions of the metric are obtained by an n-fold functional differ- 
entiation of ro[h,0,0;gni^] with respect to hfj,i,. Hereby g^i^ is kept fixed; it acts simply 
as an externally prescribed function which specifies the form of the gauge fixing condi- 
tion. Hence the functional Tq and the resulting off-shell Green functions do depend on 
(jfii^, but the on-shell Green functions, related to observable scattering amplitudes, do 
not depend on ^^j^. In this respect g^i, plays a role similar to the gauge parameter a 
in the standard approach. Remarkably, the same on-shell Green functions can be ob- 
tained by differentiating the functional r[g'ju,y] of (3.19) with respect to g/^i,, or equivalently 
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ro[/i = 0,^ = 0,^ = 0;g = g], with respect to its g argument. In this context, 'on-sheh' 
means that the metric satisfies the effective field equation Sro[g]/Sg^i, = 0. 

With (3.19) and its /c-dependent counterpart 

ffeb/ii/] = ^k[9ixu,gfiu, 0,0] (3.20) 

we succeeded in constructing a difFeomorphism invariant generating functional for gravity: 
thanks to (3.17) r[5^,^] and Tkidnu] are invariant under general coordinate transformations 
^gfif = ^vdnu- However, there is a price to be paid for their invariance: the simplified 
functional Tkidni^] does not satisfy an exact RG equation, basically because it contains 
insufficient information. The actual RG evolution has to be performed at the level of the 
functional Tk[g,g,^,^]. Only after the evolution one may set g = g, ^ = 0,^ = 0. As a 
result, the actual theory space of QEG, {A[g,g,^,^]}, consists of functionals of all four 
variables, ff/ziy, 5/*!/, ^/i, subject to the invariance condition (3.17). 

Taking a scale derivative of the regularized functional integral (3.3) and reexpressing 
the result in terms of Fj. one finds the following FRGE [11]: 



S<ri[S,{,f;5]=lTV 



hh 



2 



(3.21) 



Here F^^^ denotes the Hessian of F^ with respect to the dynamical fields h, ^, ^ at fixed g. 
It is a block matrix labeled by the fields = {hi^^, i,'^,^^,}'- 



"^^'^^'^"'"^ " vmm M(-)^^^,(y) • ^'-''^ 

(In the ghost sector the derivatives are understood as left derivatives.) Likewise, TZk is 
a block diagonal matrix with entries (^fe)^^'"' = K'^{nf^^[g])i^'''^ and H^^ = V2nf[g]. 
Performing the trace in the position representation it includes an integration J d'^x\/g{x) 
involving the background volume element. For any cutoff which is qualitatively similar to 
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the exponential cutoff tfie traces on the RHS of eq. (3.21) are well convergent, both in the 
IR and the UV. The interplay between the TZk in the denominator and the factor dfTZk in 
the numerator thereby ensures that the dominant contributions come from a narrow band 
of generalized momenta centered around k. Large momenta are exponentially suppressed. 

Besides the FRGE the effective average action also satisfies an exact integro-differential 
equation, which can be used to find the k oo limit of the average action: 

rifc-^oo[^, ^, g] = S[g + h] + S^[h; g] + 5gh[^, g] ■ (3.23) 

Intuitively, this limit can be understood from the observation that for A; — > oo all quantum 
fluctuation in the path integral are suppressed by an infinity mass-term. Thus, in this 
limit no fluctuations are integrated out and Tfc^oo agrees with the microscopic action S 

supplemented by the gauge fixing and ghost actions. At the level of the functional ^^[^(J, 
eq. (3.23) boils down to T^^^[g] = S[g]. However, as T\. involves derivatives with respect 
to h^j^i, (or equivalently gfj^n) at fixed g^j^i, it is clear that the evolution cannot be formulated 
entirely in terms of Fj^ alone. 

The background gauge invariance of F^, expressed in eq. (3.17), is of enormous practical 
importance. It implies that if the initial functional does not contain non-invariant terms, 
the flow will not generate such terms. Very often this reduces the number of terms to 
be retained in a reliable truncation ansatz quite considerably. Nevertheless, even if the 
initial action is simple, the RG flow will generate all sorts of local and non-local terms in 
Ffe which are consistent with the symmetries. 

4 Truncated flow equations 

Solving the FRGE (3.21) subject to the initial condition (3.23) is equivalent to (and in 
practice as difficult as) calculating the original functional integral over 7^,^. It is therefore 
important to devise efficient approximation methods. The truncation of theory space is 
the one which makes maximum use of the FRGE reformulation of the quantum field theory 



22 



problem at hand. 

As for the flow on the theory space {A[g, g, ^, ^]}, a still very general truncation consists 
of neglecting the evolution of the ghost action by making the ansatz 



^k[9,9,C,C] = ^k[9]+^k[9,9]+Sg{[g-g;g] + Sgh[9-9,C,C;9] , 



(4.1) 



where we extracted the classical Sg{ and Sgh from T^. The remaining functional depends 
on both g^iy and g^i,. It is further decomposed as Tk + Tk where is defined as in (3.20) 
and Ffe contains the deviations for g ^ g. Hence, by definition, rk\ci,g] = 0, and Fj. 
contains, in particular, quantum corrections to the gauge fixing term which vanishes for 
g = g, too. This ansatz satisfies the initial condition (3.23) if^ 



Tjfc^oo = S and ^k^oo = . 



(4.2) 



Inserting (4.1) into the exact FRGE (3.21) one obtains an evolution equation on the 
truncated space {74[(7,3]}: 



dt^k[9,9] 



-Tr 



-Tr 



{^-'r^k\s,9] + nr[9]y' 
[-M[g,-9]+nt[9]y'dtnt 



(4.3) 



This equation evolves the functional 



Ffe [g,g] = Ffc [g] + Sgi[g - g; g] + F^ [g,g]. 



(4.4) 



(2) 

Here F^ denotes the Hessian of rk[g,g] with respect to g^i, at fixed y^j/ and M. is given 
in eq. (3.11). 

The truncation ansatz (4.1) is still too general for practical calculations to be easily 
possible. The first truncation for which the RG flow has been found [11] is the "Einstein- 
Hilbert truncation" which retains in T^lg] only the terms / d'^x ^ and / df'x ^R, already 



See [77] for a detailed discussion of the relation between S and Sbare- 
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present in the in the classical action, with fc-dependent coupling constants, and includes 
only the wave function renormalization in r^: 

Tk[g,g] = 2K^ZNk j d^x ^ {-R + 2\k] + ^ j (fx ^r"" F^F, . (4.5) 

In this case the truncation subspace is 2-dimensional. The ansatz (4.5) contains two free 
functions of the scale, the running cosmological constant and Zjvfc or, equivalently, the 
running Newton constant Gk = G/Zi\ik- Here (5 is a fixed constant, and k = (327r(5)~"'^/^. 
As for the gauge fixing term, is given by eq. (3.9) with h^i, = g^^ — g^u replacing h^i,; 
it vanishes for g = g. The ansatz (4.5) has the general structure of (4.1) with 

Tfe = {ZNk - 1)% . (4.6) 

Within the Einstein-Hilbert approximation the gauge fixing parameter a is kept constant. 
Here we shall set a = 1 and comment on generalizations later on. 

Upon inserting the ansatz (4.5) into the flow equation (4.3) it boils down to a system 

of two ordinary differential equations for Z^vfc and X/.. Their derivation is rather technical, 

so we shall focus on the conceptual aspects here. In order to find dtZj^k and dtXk it is 

sufficient to consider (4.3) for 5^,^ = In this case the LHS of the flow equation becomes 

2k^ Jd'^x^[—RdtZ]sfk + 2dt{ZNk\k)]. The RHS is assumed to admit an expansion in 

terms of invariants Pi [gjxu] ■ In the Einstein-Hilbert truncation only two of them, J d'^x -^/g 

and Jd'^x ^R, need to be retained. They can be extracted from the traces in (4.3) by 

standard derivative expansion techniques. Equating the result to the LHS and comparing 

the coefficients of f d'^x^ and f d'^x^R, a pair of coupled differential equations for Zjvfc 

and Afc arises. It is important to note that, on the RHS, we may set g^u = g^^ only after 

(2) 

the functional derivatives of have been obtained since they must be taken at fixed 
gilv- 

As demonstrated explicitly in [27,36], this calculation can be performed without ever 
considering any specific metric g^y = g^^. This reflects the fact that the approach is 
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background covariant. The RG flow is universal in the sense that it does not depend 
on any specific metric. In this respect gravity is not different from the more traditional 
applications of the renormalization group: the RG flow in the Ising universality class, say, 
has nothing to do with any specific spin configuration, it rather reflects the statistical 
properties of very many such configurations. 

While there is no conceptual necessity to fix the background metric, it nevertheless 
is sometimes advantageous from a computational point of view to pick a specific class 
of backgrounds. Leaving ^^^^ completely general, the calculation of the functional traces 
is very hard work usually. In principle there exist well known derivative expansion and 
heat kernel techniques which could be used for this purpose, but their application is an 

(2) 

extremely lengthy and tedious task usually. Moreover, typically the operators F^. and 
TZk are of a complicated non-standard type so that no efficient use of the tabulated Seeley 
coefficients can be made. However, often calculations of this type simplify if one can 
assume that gfj,i, = g^i, has specific properties. Since the beta functions are background 
independent we may therefore restrict g^, to lie in a conveniently chosen class of geometries 
which is still general enough to disentangle the invariants retained and at the same time 
simplifies the calculation. 

For the Einstein-Hilbert truncation the most efficient choice is a family of d-spheres 
S'^{r), labeled by their radius r. These maximally symmetric backgrounds satisfy 

and, in particular, DaR^vpa = 0, so they give a vanishing value to all invariants constructed 
from g = g containing covariant derivatives acting on curvature tensors. What remains 
(among the local invariants) are terms of the form f ^P{R), where P is a polynomial 
in the Ricci scalar. Up to linear order in R the two invariants relevant for the Einstein- 
Hilbert truncation are discriminated by the S'^ metrics as the latter scale differently with 
the radius of the sphere: J^/g ^ r'^, J ^/gR ^ r'^~'^. Thus, in order to compute the beta 
functions of and Z^k it is sufficient to insert an S"^ metric with arbitrary r and to 
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compare the coefficients of r'^ and r*^"^. If one wants to do better and include the three 
quadratic invariants jR^upaR'^'^'"^ , jRfxvR^^, and /-R^, the family S'^{r) is not general 
enough to separate them; all scale like r*^""^ with the radius. 

Under the trace we need the operator [h; g\. It is most easily calculated by Taylor 
expanding the truncation ansatz, Ty;[g + h,g] = Tf;[g, g] + 0(h) + r^'^'^[^; g] + 0(h^), and 
stripping off the two ^'s from the quadratic term, 
on S'^{r) one obtains 



'i''^ = ^Jhr'^h. For g^^ the metric 



^T'^ik g] = ^ZNk^^^ J d''x I V [-D^ - 2Afe + CtR] h^'" 



with Ct = {d{d — 3) + 4)/{d{d — 1)), Cs = {d — 4:)/d. In order to partially diagonalize 
this quadratic form h^i^ has been decomposed into a traceless part h^i, and the trace 
part proportional to h^iy = h^i, + d''^g^v<i>-, g'^^h^i, = 0. Further, = g^^^D^j^Di, is the 
covariant Laplace operator corresponding to the background geometry, and R = l)/r^ 
is the numerical value of the curvature scalar on S'^{r). 

At this point we can fix the constants Zk which appear in the cutoff operators T^l'^^^ 
and TZf^ of (3.13). They should be adjusted in such a way that for every low-momentum 
mode the cutoff combines with the kinetic term of this mode to —D^ + k"^ times a constant. 
Looking at (4.8) we see that the respective kinetic terms for /i^j, and differ by a factor 
of — ((i — 2)/2d. This suggests the following choice: 



Z 



Nk • 



(4.9) 



Here {P^)^v^'^ = d g^iug^" is the projector on the trace part of the metric. For the 



traceless tensor (4.9) gives 2^^^^ 



Zj\fli.l, and for (p the different relative normalization is 



taken into account. (See ref. [11] for a detailed discussion of the subtleties related to this 
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choice.) Thus we obtain in the h and the <;!^sector, respectively: 

f/^-^rf [5,5]+7^^^)_ = Zj,A-D'' + k'R^''\-D''/k')-2\k + CTB\, (4.10) 

\ / hh \- J 

{^-'Tt\9M+nr)^^ = -^Z^,[-D' + k'R('^\-D'/k')-2-X, + CsR 

Prom now on we may set g = g and for simpUcity we have omitted the bars from the 
metric and the curvature. Since we did not take into account any renormahzation effects 
in the ghost action we set Z^^ = 1 in TZf^ and obtain 

-M + nf = -D^ + k^R^^^ {-D'^/k^) + CvR , (4.11) 

with Cv = —1/d. At this point the operator under the first trace on the RHS of (4.3) 
has become block diagonal, with the hh and (fxj) blocks given by (4.10). Both block 
operators are expressible in terms of the Laplaeian D^, in the former case acting on 
traceless symmetric tensor fields, in the latter on scalars. The second trace in (4.3) stems 
from the ghosts; it contains (4.11) with acting on vector fields. 

It is now a matter of straightforward algebra to compute the first two terms in the 
derivative expansion of those traces, proportional to J d'^x^ ~ r'' and / d'^x^R ^ r*^"^. 
Considering the trace of an arbitrary function of the Laplaeian, W{—D'^), the expansion 
up to second order derivatives of the metric is given by 

Tr[W{-D^)] = {47r)-^/hv{I)i^Qa/2[W]Jd^x^ 

+lQd/2-i[W] j d'x^R + OiR^)^ . (4.12) 

The QnS are defined as 

Qn[W\ = — - / dzz^-'Wiz) , (4.13) 
i l^'J Jo 

for n > 0, and Qo[l^] = W{0) for n = 0. The trace tr(/) counts the number of independent 
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field components. It equals 1, d, and (d— l){d + 2)/2, for scalars, vectors, and symmetric 
traceless tensors, respectively. The expansion (4.12) is easily derived using standard heat 
kernel and Mellin transform techniques [11]. 

Using (4.12) it is easy to calculate the traces in (4.3) and to obtain the RG equations 
in the form dtZjsfk = ■ ■ ■ and dt{Zj^k\) = ■ ■ ■ . We shall not display them here since it 
is more convenient to rewrite them in terms of the dimensionless running cosmological 
constant and Newton constant, respectively: 



(4.14) 



In terms of the dimensionless couplings g and A the RG equations become a system of 
autonomous differential equations 



dt9k = Pgidk, Afc) , dtXk = Pxigk, Afc) 



(4.15) 



where 



Pxig, A) = ivN - 2)A + i (47rj- " g 

X 2d{d + l)$i/2(-2A) - 8d$^/2(0) - d{d + l)777v4>^/2(-2A) 
Pg{g, X)={d-2 + riN)g. 

Here the anomalous dimension of Newton's constant rjN is given by 



l-gB2{\) 



with the following functions of the dimensionless cosmological constant: 
Si (A) ^ \ (47r)i-'^/2 ^ i)$i/^_^(-2A) - U{d - 1)^%^{-2X) 

-4d$^/2_i(0)-24$2/2(0)], 
i?2(A) ^ - i(47r)i-'^/2 \d{d + l)|.i/2_i(-2A) - 6d(d - l)4>^/2(-2A) 



(4.16) 



(4.17) 



(4.18) 



28 



The system (4.15) constitutes an approximation to a 2-dimensional projection of the RG 
flow. Its properties, and in particular the domain of apphcabihty and rehability will be 
discussed in section 5. 

The "threshold functions" $ and $ appearing in (4.16) and (4.18) are certain integrals 
involving the normalized cutoff function R^^^: 



1 ^_^ R^^\z)- zR'^^^'iz) 

Tin) Jo [z + R^^){z)+w]P ' 



1 



r(n) Jo 



dz z 



n-l 



i?(o)(z) 



[z + R^^){z) + w]P ■ 



(4.19) 



They are defined for positive integers p, and n > 0. While there are (few) aspects of 
the truncated RG flow which are independent of the cutoff scheme, i.e., independent of 
the function R^^\ the explicit solution of the flow equation requires a specific choice of 
this function. In the literature various forms of i?'^°^'s have been employed. E.g., the 
non-diff'erentiable "optimized cutofi'" [78] with R'^^\w) = (1 — w)9{l — w) allows for an 
analytic evaluation of the integrals [23] 



r(n + l) l + w; ' 



r{n + 2) 1 + w ' 



(4.20) 



Easy to handle, but disadvantageous for high precision calculations is the sharp cutoff [14] 
defined by Tlk{p^) = lim^_^^R9{l — p^/k"^), where the limit is to be taken after the 
integration. This cutoff also allows for an evaluation of the $ and $ integrals in closed 
form. Taking c/ = 4 as an example, eqs. (4.15) boil down to the following simple system 
of equations:^ 



dtXk=-{2 - m)h - — [51n(l - 2Afc) - 2C(3) + 



9t9k={^ + m) 9k , 



VN-- 



6vr + 5 fi-fc 



18 



l-2Aj 



+ 51n(l-2Afe)-C(2) + 6 



(4.21a) 
(4.21b) 
(4.21c) 



^To be precise, (4.21) corresponds to the sharp cutoff with s = 1, see [14]. 
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In order to check the scheme (in)dependence of the results it is desirable to perform 
the calculation for a whole class of /i'^^^'s. For this purpose the following one parameter 
family of exponential cutoffs has been used [13, 15, 18]: 

^^^-;^) = i^- (4-22) 

The precise form of the cutoff is controlled by the "shape parameter" s. For s = 1, 
(4.22) coincides with the standard exponential cutoff. The exponential cutoffs are suitable 
for precision calculations, but the price to be paid is that their $ and <I> integrals can be 
evaluated only numerically. The same is true for a one-parameter family of shape functions 
with compact support which was used in [13, 15]. 

Above we illustrated the general ideas and constructions underlying gravitational RG 
flows by means of the simplest example, the Einstein-Hilbert truncation. In the literature 
various extensions have been investigated. The derivation and analysis of these more 
general flow equations, corresponding to higher dimensional truncation subspaces, is an 
extremely complex and computationally demanding problem in general. For this reason 
we cannot go into the technical details here and just mention some further developments. 

(1) The natural next step beyond the Einstein-Hilbert truncation consists in generalizing 
the functional r,fc[gi], while keeping the gauge fixing and ghost sector classical, as in (4.1). 
During the RG evolution the flow generates all possible diffeomorphism invariant terms in 
Tk[g] which one can construct from g^,^. Both local and non-local terms are induced. The 
local invariants contain strings of curvature tensors and covariant derivatives acting upon 
them, with any number of tensors and derivatives, and of all possible index structures. 
The first truncation of this class which has been worked out completely [15, 16] is the 
ii^-truncation defined by (4.1) with the same Ffe as before, and the (curvature)^ action 

ffc[5]= Jd''x^[{167rGkr^[-R + 2\k]+PkR^}. (4.23) 

In this case the truncated theory space is 3-dimensional. Its natural (dimensionless) 
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coordinates are {g, A,/3), where Pk = k^^'^Pk-, and g and A defined in (4.14). Even though 
(4.23) contains only one additional invariant, the derivation of the corresponding RG 
equations is far more complicated than in the Einstein-Hilbert case. We shall summarize 
the results obtained with (4.23) [15, 16] in section 5.2. 

(2) The natural extension of the ii^-truncation consists of including all gravitational 
four-derivative terms in the truncation subspace 



k[g\ = f d^x^ [(167rGfe)-i [-R + 2h] - + ^C^^p.C^^'"' + 



(4.24) 



and adapting the classical gauge-fixing and ghost sectors to higher-derivative gravity. Here 
C^upaC^"'^'' denotes the square of the Weyl tensor, and E = Cf^^paC^'^P" - 2R^^W + 
is the integrand of the (topological) Gauss-Bonnet term in four dimensions. Using the 
FRGE (3.21), the one-loop beta functions for higher-derivative gravity have recently been 
recovered in [22,79,80], while first non-perturbative results have been obtained in [32,33]. 
The key ingredient in the non-perturbative works was the generalization of the background 
metric g from the spherical symmetric background employed in the i?^ -truncation to a 
generic Einstein background. The former case results in a flow equation for Pk = — ^ + 
while working with a generic Einstein background metric allows to find the non- 
perturbative beta functions for two independent combinations of the coupling constants 

Pk = ^ ^ ' 7fc = 7^ + — • (4.25) 

6ak oak ^CTk Ok 

The results obtained from this so-called i?^ -I- C^-truncation are detailed in section 5.3. 

(3) The first steps towards analyzing the RG flow of QEG in the ghost sector have 
recently be undertaken in [81-83]. In the flrst step, the classical ghost sector of the 
Einstein-Hilbert truncation (3.10) has been supplemented by a scale-dependent curvature 
ghost coupling C^k in the form rj^^^ = Ck J d'^x^/g R^^ where it was found that Cfe = 
constitutes a fixed point of the RG flow, with being associated with an UV-attractive 
eigendirection. The backreaction of a non-trivial ghost-wavefunction renormalization on 
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the flow of gk and was subsequently studied in [82,83], where it was estabhshed that 
the ghost-propagator is UV-suppressed by a negative anomalous dimension. Moreover, 
the phase-diagram of the ghost-improved Einstein-Hilbert truncation is almost identical 
to the one obtained without ghost-improvements shown in fig. 4. 

(4) There are also partial results concerning the gauge fixing term. Even if one makes 
the ansatz (4.5) for rfe[gi,^] in which the gauge fixing term has the classical (or more 
appropriately, bare) structure one should treat its prefactor as a running coupling: a = a^. 
The beta function of a has not been determined yet from the FRGE, but there is a simple 
argument which allows us to bypass this calculation. 

In non-perturbative Yang-Mills theory and in perturbative quantum gravity a = ak = 
is known to be a fixed point for the a evolution. The following reasoning suggests that 
the same is true within the non-perturbative FRGE approach to gravity. In the standard 
functional integral the limit a — > corresponds to a sharp implementation of the gauge 
fixing condition, i.e., cxp(— S'gf) becomes proportional to The domain of the / Vh^^, 

integration consists of those h^^'s which satisfy the gauge fixing condition exactly, = 0. 
Adding the IR cutoff at k amounts to suppressing some of the /i^j^ modes while retaining 
the others. But since all of them satisfy = 0, a variation of k cannot change the domain 
of the hfj,,^ integration. The delta functional S[Ffj] continues to be present for any value of 
k if it was there originally. As a consequence, a vanishes for all k, i.e., a = is a fixed 
point of the a evolution [84] . 

Thus we can mimic the dynamical treatment of a running a by setting the gauge fixing 
parameter to the constant value a = 0. The calculation for a = is more complicated 
than at a = 1, but for the Einstein-Hilbert truncation the a-dependence of Pg and j3x, for 
arbitrary constant a has been found in [13,85]. The i?^ -truncations could be analyzed only 
in the simple a = 1 gauge, but the results from the Einstein-Hilbert truncation suggest 
the UV quantities of interest do not change much between a = and a = 1 [13, 15]. 

(5) In refs. [28-30] gravitational RG flows have been explored in an approximation that 
goes beyond a truncation of theory space. Here only the subsector of the basic path 
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integral over the conformal degrees of freedom has been considered while all others were 
omitted. This leads to a scalar-like theory to which the same apparatus underlying the 
analysis of the full theory has been applied (average action, background decomposition). 
Remarkably, this conformally reduced gravity (contrary to a standard 4-dimensional 
scalar theory) possesses a NGFP and a RG flow that is qualitatively similar to that of full 
QEG. It was possible to establish the existence of the NGFP on an infinite dimensional 
theory space consisting of arbitrary potentials for the conformal factor. These somewhat 
unexpected results find their explanation [28, 29] by noting that the quantization scheme 
based upon the (conformal reduction of the) gravitational average action is "background 
independent" in the sense that no special metric (flat space, etc.) plays a distinguished 
role. 

(6) Up to now we considered pure gravity. As for as the general formalism, the inclusion 
of matter fields is straightforward. The structure of the flow equation remains unaltered, 
except that now F^^^ and TZk are operators on the larger Hilbert space of both gravity and 
matter fluctuations. In practice the derivation of the projected RG equations can be 
quite a formidable task, however, the difficult part being the decoupling of the various 

(2) 

modes (diagonalization of Fj^ '^) which in most calculational schemes is necessary for the 
computation of the functional traces. Various matter systems, both interacting and non- 
interacting (apart from their interaction with gravity) have been studied in the literature 
[12,86,87]. A rather detailed analysis has been performed by Percacci et al. In [12,21] 
arbitrary multiplets of free (massless) fields with spin 0, 1/2, 1 and 3/2 were included. In 
[21] an interacting scalar theory coupled to gravity in the Einstein-Hilbert approximation 
was analyzed, and a possible solution to the triviality and the hierarchy problem [88] was 
a first application in this context. 

(7) At the perturbative one-loop level, the fiow of possibly non-local form factors 
appearing in the curvature expansion of the effective average action has been studied 
in [89,90]. For a a minimally coupled scalar field on a 2-dimensional curved space-time, the 
fiow equation for the form factor in J d'^x^Rck{A)R correctly reproduces the Polyakov 
effective action, while in d = 4 this ansatz allows to recover the low energy effective action 
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as derived in the effective field tfieory framework [91]. 

(8) As yet, almost all truncations studied arc of the "single metric" type where Tk depends 
on via the gauge fixing term only. The investigation of genuine bimetric truncations 
with a nontrivial dependence on both g^i, and g/^i, started only recently in [92,93]. Future 
work on the verification of the NGFP has to go in this direction clearly. 

(9) In ref. [94] a first investigation of Lorentzian gravity was performed in a 3+1 split 
setting. In this case the flow equation (1.1) has been formulated in terms of the ADM- 
decomposed metric degrees of freedom, which imprints a foliation structure on space-time 
and provides a preferred "time-direction" . The resulting FRGE depends on an additional 
parameter e, which encodes the signature of the space-time metric. For a fixed truncation 
the resulting Lorentzian renormalization group flow turned out almost identical to the one 
obtained in the Euclidean case. 

(10) A first step towards new gravitational field variables, different from the metric, 
was taken in [95]. Employing the vielbein and spin connection as independent variables, 
the question of Asymptotic Safety was reconsidered. In principle it is conceivable that an 
Einstein-Cartan type quantum field theory is inequivalent to QEG. However, the actual 
calculations support the conjecture that "Quantum Einstein-Cartan Gravity" is asymp- 
totically safe, too.^ 

5 Average action approach to Asymptotic Safety 

Based on the exact flow equation (1.1), we now implement the ideas of the Asymptotic 
Safety construction at the level of explicitly computable approximate RG flows on trun- 
cated theory spaces. A summary of the truncations explored to date is provided in fig. 
3. For a detailed derivation of the beta functions we refer to [11] (Einstein-Hilbert trun- 
cation), [15] (i?^-truncation), [26] (/(i?)-truncation), and [33] for the i?^ + C^-truncation, 
respectively. 

®Also see [96] for a perturbative analysis of the running Immirzi parameter. 
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Figure 3: Overview of the various truncations employed in the systematic exploration of the theory 
space of QEG. The lines indicate the interaction monomials contained in the various truncation 
ansatzc for f ^[g], cq. (4.4). All truncations have confirmed the existence of a non-trivial UV fixed 
point of the gravitational RG flow. 

5.1 The Einstein-Hilbert truncation 

The Einstein-Hilbert truncation (4.5) constitutes the most prominent truncation studied 
to date [11,13,14,18,23,27,36]. In [14] the corresponding RG equations (4.15) have been 
analyzed in detail, using both analytical and numerical methods. In particular all RG 
trajectories have been classified, and examples have been computed numerically. The 
most important classes of trajectories in the phase portrait on the ^-A— plane are shown 
in fig. 4. Notably, all cutoflFs tested to date confirm this picture at least qualitatively. 

The RG flow is found to be dominated by two fixed points {g*,\*): the GFP at 
g* = \* = 0, and a NGFP with g* > and A* > 0. There arc three classes of trajectories 
emanating from the NGFP: trajectories of Type la and Ilia run towards negative and 
positive cosmological constants, respectively, and the single trajectory of Type Ila ("sepa- 
ratrix") hits the GFP for k 0. The high momentum properties of QEG are governed by 
the NGFP; for A; — > oo, in fig. 4 all RG trajectories on the half-plane g > run into this 
point. The fact that at the NGFP the dimensionless coupling constants gk,^k approach 
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Figure 4: RG flow in the g-A— plane. The arrows point in the direction of increasing coarse 
graining, i.e., of decreasing k. (From [14].) 

constant, non-zero values then implies that the dimensionful quantities run according to 



Hence for ^ oo and d > 2 the dimensionful Newton constant vanishes while the cosmo- 
logical constant diverges. 

Thus the Einstein-Hilbert truncation does indeed predict the existence of a NGFP with 
exactly the properties needed for the Asymptotic Safety construction. Clearly the crucial 
question now is whether this NGFP is the projection of a fixed point in the exact theory or 
whether it is merely the artifact of an insufficient approximation. We now summarize the 

properties of the NGFP established within the Einstein-Hilbert truncation. All findings 
mentioned below are independent pieces of evidence pointing in the direction that QEG 
is indeed asymptotically safe in four dimensions. Except for point (5) all results refer to 



(1) Universal existence: The non-Gaussian fixed point exists for all cutoff schemes and 



Gk = 9*1^ 



(5.1) 



d = 4. 



36 



shape functions implemented to date. It seems impossible to find an admissible cutoff 
which destroys the fixed point in d = 4. This result is highly non-trivial since in higher 
dimensions {d > 5) the existence of the NGFP depends on the cutoff chosen [14]. 

(2) Positive Newton constant: While the position of the fixed point is scheme de- 
pendent, all cutoffs yield positive values of g* and A*. A negative g* might have been 
problematic for stability reasons, but there is no mechanism in the flow equation which 
would exclude it on general grounds. 

(3) Stability: For any cutoff employed the NGFP is found to be UV attractive in both 
directions of the A-^— plane. Linearizing the flow equation according to eq. (2.6) we obtain 
a pair of complex conjugate critical exponents 9i = 0| with positive real part 6' and 
imaginary parts ±6". In terms oi t = liL{k/ko) the general solution to the linearized flow 
equations reads 

{h,9kf = {X*,g*f + 2f^[ReC cos{e"t) +ImC sm{e"t)]ReV 

+ [ReC sin(^"t) - ImC cos Imyje"^'* . (5.2) 

with C = Ci = {C2)* an arbitrary complex number and V = = {V'^)* the right- 
eigenvector of B with eigenvalue —9i = — ^1- -^1- (2-6) implies that, due to the positivity 
of 0', all trajectories hit the fixed point as i is sent to infinity. The non- vanishing imaginary 
part 6" has no impact on the stability. However, it influences the shape of the trajectories 
which spiral into the fixed point for A; — > 00. Thus, the fixed point has the stability 
properties needed in the Asymptotic Safety scenario. 

Solving the full, non-linear flow equations [14] shows that the asymptotic scaling region 
where the linearization (5.2) is valid extends from k = 00 down to about k ^ mpi with 
the Planck mass deflned as mp\ = Gq^^^. Here mpi plays a role similar to Aqcd in QCD: 
it marks the lower boundary of the asymptotic scaling region. We set ko = mpi so that 
the asymptotic scaling regime extends from about t = to t = 00. 
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(4) Scheme- and gauge dependence: Analyzing the cutoff scheme dependence of 0' , 
9", and g*X* as a measure for the rehabihty of the truncation, the critical exponents were 
found to be reasonably constant within about a factor of 2. For a = 1 and a = 0, for 

instance, they assume values in the ranges 1.4 < 9' < 1.8, 2.3 < ^" < 4 and 1.7 < 9' < 2.1, 
2.5 < ^" < 5, respectively. The universality properties of the product g*X* are even more 
impressive. Despite the rather strong scheme dependence of g* and A* separately, their 
product has almost no visible s-dependence for not too small values of s. Its value is 



9 A 



0.12 for a = 1 

(5.3) 

0.14 for a = 0. 



The difference between the "physical" (fixed point) value of the gauge parameter, a = 0, 
and the technically more convenient a = 1 are at the level of about 10 to 20 percent. 

(5) Higher and lower dimensions: The beta functions implied by the FRGE are 
continuous functions of the space-time dimensionality and it is instructive to analyze 
them for d ^ 4. In rcf. [11] it has been shown that for d = 2 + e, |e| <C 1, the FRGE 
reproduces Weinberg's [9] fixed point for Newton's constant, g* = ^e, and also supplies 
a corresponding fixed point value for the cosmological constant. A* = — ^<I>^(0)e, with 
the threshold function given in (4.19). For arbitrary d and a generic cutoff the KG flow 
is quantitatively similar to the 4-dimensional one for all d smaller than a certain critical 
dimension dcrit; above which the existence or non-existence of the NGFP becomes cutoff- 
dependent. The critical dimension is scheme dependent, but for any admissible cutoff 
it lies well above d = 4. As d approaches dcrit from below, the scheme dependence of 
the universal quantities increases drastically, indicating that the -R-truncation becomes 
insufficient near dcrit- 

In fig. 5 we show the d-dependence of g*, A*, 0', and 9" for two versions of the sharp 
cutoff (with s = 1 and s = 30, respectively) and for the exponential cutoff with s = 1. For 
2 + s < d < 4 the scheme dependence of the critical exponents is rather weak; it becomes 
appreciable only near d « 6 [14]. Fig. 5 suggests that the Einstein-Hilbert truncation in 
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Figure 5: Comparison of X*,g*,9' and 9" for different cutoff functions in dependence of the 
dimension d. Two versions of the sharp cutoff (sc) and the exponential cutoff with s = 1 (Exp) 
have been employed. The upper line shows that for 2 + e < d < 4 tht' cutoff scheme dependence of 
the results is rather small. The lower diagram shows that increasing d beyond about 5 leads to a 
significant difference in the results for 6' , 9" obtained with the different cutoff schemes. (From [14].) 

d = 4 performs almost as well as near d = 2. Its validity can be extended towards larger 
dimensionalities by optimizing the shape function [23]. 

5.2 /(i?)-type truncations 

The ultimate justification of a given truncation consists in checking that if one adds further 
terms to it, its physical predictions remain robust. The first step towards testing the 
robustness of the Einstein-Hilbert truncation near the NGFP against the inclusion of 
other invariants has been taken in refs. [15, 16] where the beta functions for the three 
generalized couplings g, A and /3 entering into the i?^-truncation of eq. (4.23) have been 
derived and analyzed. 
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Subsequently, the truncated theory space has been extended to arbitrary functions 
of the Ricci scalar in [25-27, 35] . In this truncation ansatz S'gf and Sgt are taken to be 
classical while, in the language of cq. (4.1), 

ffcb] = J d^x^fkiR) , fk[9,g] = 0. (5.4) 

Substituting this ansatz into the flow equation (1.1) results in a rather complicated partial 
differential equation governing the scale-dependence of fk{R) [26]. Based on this equation, 
the search for the NGFP on truncation subspaces involving higher powers of the curvature 
scalar reduces to an algebraic problem. Substituting the ansatz^ 

N 

h{R) = Y,Un{k)k\R/k^r , iVGN, (5.5) 

n=0 

and expanding the resulting equation in powers of R allows to extract the non-perturbative 
beta functions for the dimensionless couplings Un{k), 

kdkUn{k) = PuAuo,--- ,un), n = 0, (5.6) 

The fixed point conditions /3u„(nQ, • • • , u*^) = 0,n = 0, . . . , N can then be solved numer- 
ically. Notably, the inclusion of higher-derivative terms provided crucial evidence that 
UV-critical hypersurface of the NGFP known from the Einstein-Hilbert truncation has a 
finite dimension, indicating that the Asymptotic Safety scenario is predictive. We shall 
now summarize the central results obtained within this class of truncations. 

(1) Position of the fixed point (i?^): Also with the generalized truncation (4.23) the 

NGFP is found to exist for all admissible cutoffs. Fig. 6 shows its coordinates {X* , g* , /3*) 

for the family of shape functions (4.22). For every shape parameter s, the values of 

A* and g* are almost the same as those obtained with the Einstein-Hilbert truncation. 

^The Einstein-Hilbert truncation discussed in section 4 corresponds to setting fk{R) = 
(167rGfc)^^ (— -R-|-2Afe) and using (4.6) rather than setting r[g,g] = 0. Comparing the Einstein-Hilbert 
ansatz to (5.5) shows ui{k) = — (167rgfc)~^, so that a negative wi(fc) actually corresponds to a positive 
Newton's constant. 
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(a) (b) 

Figure 6: (a) g* , A*, and g*\* as functions of s for 1 < .s < 5, and (b) (3* as a function of s for 
1 < s < 30, using the family of exponential shape functions (4.22). (From ref. [16].) 

In particular, the product g*\* is constant with a very high accuracy. For s = 1, for 
instance, one obtains (A*,(/*) = (0.348,0.272) from the Einstein-Hilbert truncation and 
{X* , g* , f3*) = (0.330,0.292,0.005) from the generalized truncation. It is quite remarkable 
that (3* is always significantly smaller than A* and g*. Within the limited precision of 
our calculation this means that in the 3-dimensional parameter space the fixed point 
practically lies on the A-^— plane with ^ = 0, i.e., on the parameter space of the pure 
Einstein-Hilbert truncation. 

(2) Eigenvalues and -vectors (i?^): The NGFP of the i2^-truncation proves to be UV 
attractive in any of the three directions of the (A, ^f, /?)— space for all cutoffs used. The 
linearized flow in its vicinity is always governed by a pair of complex conjugate critical 
exponents 9i = 9' + \6" = 62 with ^' > and a single real, positive critical exponent 
9z > 0. For the exponential shape function with s = 1, for instance, we find 9' = 2.15, 
9" = 3.79, 93 = 28.8. The first two are again of the spiral type while the third one is a 
straight line. 

For any cutoff, the numerical results have several quite remarkable properties. They 
all indicate that, close to the NGFP, the RG flow is rather well approximated by the pure 
Einstein-Hilbert truncation. 

(a) The eigenvectors associated with the spiraling directions span a plane which virtually 
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Figure 7: Trajectory of the linearized flow equation obtained from the i?^-truncation for 1 < 
t = ln(fc/fco) < oo. In (b) we depict the eigendirections and the "box" to which the trajectory is 
confined. (Prom ref. [16].) 

coincides with the ^-A— subspace at /3 = 0, i.e., with the parameter space of the Einstein- 
Hilbert truncation. As a consequence, the corresponding normal modes are essentially the 
same trajectories as the "old" normal modes already found without the i?^-term. Also 
the corresponding 0'- and ^"-values coincide within the scheme dependence. 

(b) The new eigenvalue 63 introduced by the i?^-term is significantly larger than 0'. When 
a trajectory approaches the fixed point from below {t 00), the "old" normal modes are 
proportional to exp{—9't), but the new one is proportional to exp(— ^3*), so that it decays 
much quicker. For every trajectory running into the fixed point we find therefore that 
once t is sufficiently large the trajectory lies entirely in the /3 = 0-plane practically. Due 
to the large value of ^3, the new scaling field is very ''relevant". However, when we start 
at the fixed point (t = 00) and lower t it is only at the low energy scale k « mp\ (t « 0) 
that exp(— 03i) reaches unity, and only then, i.e., far away from the fixed point, the new 
scaling field starts growing rapidly. 

Thus very close to the fixed point the RG flow seems to be essentially 2-dimensional, 
and that this 2-dimensional flow is well approximated by the RG equations of the Einstein- 
Hilbert truncation. In fig. 7 we show a typical trajectory which has all three normal modes 
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Figure 8: (a) 6' = ReOi and 6" = Im^i, and (b) ^3 as functions of s, using tlie family of 
exponential shape functions (4.22). (From [15].) 

excited with equal strength. All the way down from A; = 00 to about k = mpi it is confined 
to a very thin box surrounding the (3 = 0-plane. 

(3) Scheme dependence (i?^): The scheme dependence of the critical exponents and 
of the product g*X* turns out to be of the same order of magnitude as in the case of the 
Einstein-Hilbcrt truncation. Fig. 8 shows the cutoff dependence of the critical exponents, 
using the family of shape functions (4.22). For the cutoffs employed 9' and 9" assume 
values in the ranges 2.1 < 6' < 3.4 and 3.1 < 9" < 4.3, respectively. While the scheme 
dependence of 9" is weaker than in the case of the Einstein-Hilbert truncation one finds 
that it is slightly larger for 9'. The exponent 9s suffers from relatively strong variations 
as the cutoff is changed, 8.4 < ^3 < 28.8, but it is always significantly larger than 9'. The 
product g*X* again exhibits an extremely weak scheme dependence. Fig. 6(a) displays 
g*X* as a function of s. It is impressive to see how the cutoff dependences of g* and A* 
cancel almost perfectly. Fig. 6(a) suggests the universal value g*X* ~ 0.14. Comparing 
this value to those obtained from the Einstein-Hilbert truncation we find that it differs 
slightly from the one based upon the same gauge a = 1. The deviation is of the same 
size as the difference between the a = 0- and the a = 1-results of the Einstein-Hilbert 
truncation. 

(4) Dimensionality of 5uv (R^)'- According to the canonical dimensional analysis, the 
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0.955£ 
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Table 1: Fixed point coordinates and critical exponents of the ii^- truncation in 2 + e dimensions. 
The negative value ^3 < implies that <Suv is an (only!) 2-dimensional surface in the 3-dimensional 
theory space. 

(curvature) "-invariants in 4 dimensions are classically marginal forn = 2 and irrelevant for 
n > 2. The results for 63 indicate that there are large non-classical contributions so that 
there might be relevant operators perhaps even beyond n = 2. With the i?^— truncation 
it is clearly not possible to determine their number Auv in d = 4. However, as it is 
hardly conceivable that the quantum effects change the signs of arbitrarily large (negative) 
classical scaling dimensions, Auv should be finite [9]. 

The first confirmation of this pictTirc came from the /^'^-calculation in d = 2 -|- e where 
the dimensional count is shifted by two units. In this case we find indeed that the third 
scaling field is irrelevant for any cutoff employed, ^3 < 0. Using the £-expansion the 
corresponding numerical results for selected values of the shape parameter s are presented 
in table 1. 

For all cutoffs used we obtain three real critical exponents, the first two are positive 
and the third is negative. This suggests that in d = 2 -|- e the dimensionality of Sjjy could 
be as small as Auv = 2 and characterized by only two free parameters, the renormalized 
Newton constant Go and the renormalized cosmological constant Aq, for instance. 

(5) Position of the non-Gaussian fixed point {f{R),d = 4): The beta functions 
(5.6) also give rise to a NGFP with g* > 0, A* > 0, whose A^-dependent position is shown 
in table 2. In particular the product g*X* = uq/{32tt{ui)'^) displayed in the last column 
is remarkably constant. It is in excellent agreement with the Einstein-Hilbcrt truncation 
(5.3), and the i?^-truncation, fig. 6(a). Only the value obtained in the case N = 2 shows 
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Table 2: Location of the NGFP obtained within the /(-R)-truncation by expanding the partial 
differential equation in a power series in R up to order , including fc-dependent dimensionless 
coupling constants Un{k), n = 0, • • • , A'^. (From [26].) 
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Table 3: Stability coefficients of the NGFP for increasing dimension + 1 of the truncation 
subspace. The first two critical exponents are a complex pair 9± = 9' ± W . (From [26].) 

a mild deviation from the i?^— computations, which can, most probably, be attributed to 
the use of a different gauge-fixing procedure, cutoff shape function, and ansatz for V^^. 

(6) Dimensionality of S\jy {f{R),d = 4): The critical exponents resulting from the 
stability analysis of the NGFP emerging from the /(i?)-truncatlon are summarized in table 
3. In particular we find that only three of the eigendirections associated to the NGFP are 
relevant, i.e., UV attractive. Including higher derivative terms i?", n > 3 in the truncation 
creates irrelevant directions only. Thus, as In the R^- and + C^-truncation, Auv = 3, 
i.e., the UV critical surface associated to the fixed point Is a 3-dimensional submanifold 
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in the truncated theory space. Its dimensionahty is stable with respect to increasing TV. 
The RG trajectories tracing out this surface are determined by fixing the three relevant 
couplings. These describe in which direction tangent to 5uv the trajectory flows away 
from the fixed point. All remaining couplings, the irrelevant ones, are predictions from 
Asymptotic Safety. In [27] these results have been extended up to = 8, providing 
even stronger evidence for the robustness of the RG flow under the inclusion of further 
invariants. 

5.3 The + C^-truncation in d = 4 

The extension of the i?^-truncation by including non-scalar curvature terms in the trun- 
cation subspace has been investigated in [32,33]. This setup is based on the ansatz (4.24), 
and includes the characteristic features of higher-derivative gravity as, e.g., a fourth-order 
propagator for the helicity 2 states. Moreover, it provides an important spring-board for 
understanding the role of the counterterms arising in the perturbative quantization of 
gravity in the Asymptotic Safety program. 

Including tensor structures like C^j,i,p„C^^^f"^ in the truncation subspace requires the 
generalization of the background metrics g. As discussed at the end of section 4, choosing 
the class of maximally symmetric metrics on 5*^ as a background considerably simplifies 
the evaluation of the truncated flow equation, but comes with the drawback that the flow 
is projected on interaction terms built from the Ricci scalar only. Thus the inclusion of the 
remaining four-derivative operators requires the use of a more general class of backgrounds. 
Ideally, this new class is generic enough to disentangle the coefficients multiplying B? and 
the tensorial terms, and, most importantly, simple enough to avoid the appearance of 
non-minimal higher-derivative differential operators inside the trace. While the maximally 
symmetric backgrounds used up to now are insufficient in the former respect, a generic 
compact Einstein background (without Killing or conformal Killing vectors and without 
boundary for simplicity), satisfying R^y = g^j_u, is sufficient to meet both criteria and 
allows one to determine the non-perturbative beta functions of the linear combinations 
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(4.25). 

Surprisingly, projecting the flow equation resulting from the ansatz (4.24) onto an 
generic Einstein background, the differential operators appearing on its RHS organize 
themselves into second order differential operators of the Lichnerowicz form 

= -5'<^M. - 2i?^V'/'a/3 > AiLcl)^, = -D^^-R^,r , AoL<^ = -L'20, (5.7) 

which commute with all the other curvature terms inside the trace. This feature makes the 
traces amenable to standard heat kernel techniques for minimal second order differential 
operators. The resulting beta functions for the dimensionless coupling constants \k,9k,Pk, 
and 7fc are somewhat involved, so that we only highlight their main properties 

(1) Existence of the NGFP: The beta functions of the + C^-truncation also give 
rise to a NGFP with positive Newtons and cosmological constant 

g* = 1.960 , A* = 0.218 , P* = 0.008 , 7* = -0.005 , g*X* = 0.427 . (5.8) 

The finite values for P* and 7* also imply a finite value of a*, via eq. (4.25). This should 
be contrasted to the one-loop result a* = obtained from the perturbative quantization of 
fourth-order gravity. Thus the non-perturbative corrections captured by the FRGE shift 
the fixed point underlying the asymptotic freedom obtained within perturbation theory 
to the NGFP featuring in the Asymptotic Safety program. 

(2) Stability properties: An important characteristics of the NGFP are its stability 
properties. Linearizing the RG flow at the fixed point (5.8) along the lines of section 2, 
the stability coefficients are found as 

^0 = 2.51, 01 = 1.69, 02 = 8.40, 63 = -2.11. (5.9) 

We observe that the inclusion of the C^-coupling leads to real stability coefficients. This is 
in contrast to the complex stability coefficients and the corresponding spiraling approach 
of the RG flow characteristic for /(i?)-type truncations. 
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Moreover, eq. (5.9) further indicates that the "new" direction added by extending the 
i?^-truncation is UV-repulsive. Thus, in this 4-dimensional truncation space, Sjjy remains 
3-dimensional, as in the case of the i?^-truncation. The condition for a trajectory being 
inside the UV-critical surface then imposes one constraint between the coupUng constants. 
In the hnear regime at the NGFP, this constraint can be used to express 7^, say, in terms 
of the other couphng constants contained in the ansatz 

7fe = -0.116 + 0.030 Afe g^^ + 0.049 5"^ + 11.06^^ . (5.10) 

Comparing the results (5.8) and (5.9) to their counterparts in the Einstein-Hilbert 
and i?^-truncation summarized in the previous subsections, we conclude that the C^-term 
leads to a moderate shift on fixed point structure. In particular the universal coupling 
g*X* turns out to be enhanced by a factor of three, while the stability coefficients of the 
fixed point turn out to be real. Thus the C^-term significantly infiuences the RG flow 
of the theory. It effects the fixed point structure more drastically then the inclusion of 
the i?^-term or working with different cutoff-schemes or gauge-fixing functions within the 
Einstein-Hilbert action. 

(3) The role of perturbative counterterms: An interesting twist arises from supple- 
menting the truncation ansatz (4.24) with a free scalar field. From perturbative viewpoint 
the two-derivative terms of the resulting action constitute the prototypical example of a 
gravitational theory which is perturbatively non-renormalizable at one loop, with the four- 
derivative terms being the corresponding non-renormalizable on-shell counterterms [2]. 
At the same time, it has been shown in [21] that the RG flow obtained from the two- 
derivative truncation possesses a NGFP whose properties are very close to the one found 
in the Einstein-Hilbert truncation of pure gravity. Thus this setup provides a valuable 
laboratory, where the influence of perturbative counterterms on the NGFP underlying the 
Asymptotic Safety scenario can be studied explicitly. 

The beta function for the "i?^ + + scalar" -truncation can be computed completely 
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analogous to the case of pure gravity. Notably, they also give rise to an NGFP [33] 

5* = 2.279, A* =0.251, /3* = 0.010, 7* = -0.0043, /A* = 0.571 , (5.11) 

with stability coefficients 

00 = 2.67, 01 = 1.39, 02 = 7.86, 03 = -1-50. (5.12) 

Comparing these to the pure gravity results (5.8), (5.9), we observe that the gravity- 
matter fixed point has strikingly similar properties. The universal couplings and stability 
coefficients change their values by roughly 25%. Thus the beta functions are dominated 
by the contribution from the gravitational sector. This result lends strong support to the 
assertion that perturbative counterterms do not play a distinguished role in the Asymptotic 
Safety program. In particular they are not fatal for the NGFP, so that the gravity-matter 
theory remains asymptotically safe despite the inclusion of an perturbative conterterms in 
the truncation subspace. 

Summary: The above results strongly suggest that the non-Gaussian fixed point occur- 
ring in the Einstein-Hilbert truncation is not a truncation artifact but rather the projection 
of a fixed point in the exact theory space. The fixed point and all its qualitative properties 
are stable against variations of the cutoff and the inclusion of further invariants in the trun- 
cation. It is particularly remarkable that within the scheme dependence the additional 
i2^-term has essentially no impact on the fixed point. Moreover, truncations involving 
higher-order polynomials in R or the tensor structure fully confirm the general picture 
suggested by the simple Einstein-Hilbert truncation and provide a strong indication that 
the corresponding quantum field theories are characterized by a finite number of free pa- 
rameters. We interpret the above results and their mutual consistency as quite non-trivial 
indications supporting the conjecture that 4-dimensional QEG indeed possesses a RG fixed 
point with precisely the properties needed for Asymptotic Safety. 
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6 The multifractal properties of QEG space-times 



An intriguing consequence arising from the scale-dependence of the gravitational effective 
action is that the QEG space-times at short distances develop fractal properties [13, 15, 
49,60]. These manifest themselves at the level of the graviton propagator at high energies 
or in the diffusion processes of test particles. One of the striking conclusions reached in 
refs. [13, 15] was that the effective dimensionality of space-time equals 4 at macroscopic 
distances {£ ':$> ipi) but, near £ £pi, it gets dynamically reduced to the value 2. The 
remainder of this review is dedicated to the discussion of the arguments that underlying 
this conclusion. 



6.1 The origin of fractality: scale-dependent metrics 

As we have seen, the effective average action rfe[5^i^] defines a continuous family of effective 
field theories, where all quantum fluctuations with momenta larger than k have been 
integrated out. Intuitively speaking, the solution {g^v)^ of the scale dependent field 
equation 

k5>J=0 (6-1) 



can be interpreted as the metric averaged over (Euclidean) space-time volumes of a linear 
extension £ which typically is of the order of 1/k. Knowing the scale dependence of r^, 
i.e., the renormalization group trajectory k i-> F^, we can in principle follow the solution 
{9nv)k from the ultraviolet {k oo) to the infrared {k 0). 

It is an important feature of this approach that the infinitely many equations of (6.1), 
one for each scale k, are valid simultaneously. They all refer to the same physical system, 
the "quantum space-time", but describe its effective metric structure on different scales. 
An observer using a "microscope" with a resolution k~^ will perceive the universe to be 
a Riemannian manifold with metric {g^u)/^- At every fixed k, (5^1/)^ is a smooth classical 
metric. But since the quantum space-time is characterized by the infinity of equations 
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(6.1) with k = 0, - ■ ■ , oo it can acquire very nonclassical and in particular fractal features. 

Let us describe more precisely what it means to "average" over Euclidean space-time 
volumes. The quantity we can freely tune is the IR cutoff scale k. The "resolving power" 
of the microscope, henceforth denoted £, is in general a complicated function of k. (In flat 
space, £ ^ 1/k.) In order to understand the relationship between £ and k we must recall 
some steps from the construction of Tklgui/] in section 3. 

The IR cutoff of the average action is implemented by first expressing the functional 
integral over all metrics in terms of eigenmodes of D^, the covariant Laplacian formed with 
the aid of the background metric 5^;^. Then a suppression term is introduced which damps 
the contribution of all — 5^-modes with eigenvalues smaller than k^. Following the steps 
of section 3 this leads to the scale dependent functional ^k[9ij,i/',9ij,i/], and the action with 
one argument is again obtained by equating the two metrics: rjfc[5^i^] = T^lgi^i^; g^^ = g^i^]. 
This is this action which appears in (6.1). Because of the identification of the two metrics 
we see that it is basically the eigenmodes of D"^ = D"^, constructed from the argument of 
rfcis]) which are cut off at k'^. Since {gf^u)j^ is the corresponding stationary point, we can 
say that the metric (5^,/)^ applies to the situation where only the quantum fluctuations 
of —D'^{(^g^v)^) with eigenvalues larger than k"^ are integrated out. Therefore there is a 
complicated interrelation between the metric and the scale at which it provides an effective 
description: The covariant Laplacian which ultimately decides about which modes are 
integrated out is constructed from the "on shell" configuration {g^v)j^-, so it is /c-dependent 
by itself already. 

Prom these remarks it is clear now how to obtain the "resolving power" £ for a given 
A;, at least in principle. We start from a fixed RG trajectory k i-> Fjt, derive its ef- 
fective field equations at each k, and solve them. The resulting quantum mechanical 
counterpart of a classical space-time is specified by the infinity of Riemannian metrics 
{{gij,i')^\h = 0, • ■ • ,00}. While the totality of these metrics contains all physical infor- 
mation, the parameter k is only a book keeping device a priori. In a second step, it 
can be given a physical interpretation by relating it to the (proper) length scale of the 
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averaging procedure: One constructs the Laplacian —D'^{(^g^v')^, diagonalizes it, looks 
how rapidly its A;^-eigenfunction varies, and "measures" the length I of typical variations 
with the metric {gnv)j^ itself. By solving the resulting I = l{k) for k = k{£) we can in 
principle reinterpret the metric {gfi,/)f^ as referring to a microscope with a known position 
and direction dependent resolving power. The price we have to pay for the background 
independence is that we cannot freely choose i directly but rather k only. 

We now illustrate this procedure at the level of the Einstein-Hilbert truncation dis- 
cussed in section 5.1. Without matter, the corresponding field equations happen to co- 
incide with the ordinary Einstein equation, but with Gk and Afc replacing the classical 
constants 

R^A{9)k) = ^>^k{9,.),. (6.2) 

It is easy to make the /c-dcpcndcncc of ((7^;/)^ explicit. Picking an arbitrary refer- 
ence scale ko we may rewrite (6.2) as [Xko/Xk] R'^u{{g)f^) = ^r^^ko^ii- Since R^u{cg) = 

R^v{g) for any constant c > 0, this relation implies that the average metric and its 
inverse scale as 

{9A^))k = [K/h]{g^^uix))ko , {g'"'ix))k = [h/K]{g''''ix))ko ■ (6.3) 

These relations are valid provided the family of solutions considered exists for all scales 
between ko and k, and Afc has the same sign always. 

Denoting the Laplace operators corresponding to the metrics {gny)k and {gnu)ko by 
A(fe) and A(fco), respectively, these relations imply 

Aik) = [Xk/ho] A(A:o) . (6.4) 

At this stage, the following remark is in order. In the asymptotic scaling regime 
associated with the NGFP the scale-dependence of the couplings is determined by the 
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fixed point, see eq. (5.1). Choosing Afeg in the classical regime, this imphes in particular 

{g^u{x))k oc k~'^ {k^oo). (6.5) 

This asymptotic relation is actually an exact consequence of Asymptotic Safety, which 
solely relies on the scale-independence of the theory at the fixed point. This can be seen 
as follows. The complete effective average action has the structure (2.2). If Ua{k) has the 
canonical dimension da the corresponding dimensionless couplings are Ua{k) = k~'^°' Ua{k) 
and we have 

a a 

In the second equality we used that Pa [c^ g^v] = c'^" Pa [dtiv] for any c > since Pa has 
dimension —da- If the theory is asymptotically safe at the exact level, all Ua{k) approach 
constant values u* for A; ^ oo: 

^k-^oo[gfiv\ = ^U*aPa[k^ Qfiu]- (6.7) 
a 

Obviously this functional depends on /c^ and g^^^ only via the combination k"^ g^^. There- 
fore the solutions of the corresponding field equation, (5^^)^, scale proportional to 
and this is exactly the scaling behavior (6.5). 

The fact that (6.3) implies a fractal structure of space-time at short distances can be 
argued as follows. Since in absence of dimensionful constants of integration is the only 
quantity in this equation which sets a scale, every solution to (6.2) has a typical radius 
of curvature rc{k) oc l/-\/Afc. (For instance, the maximally symmetric S'^-solution has the 
radius rg = r = y3/Xk.) If we want to explore the space-time structure at a fixed length 
scale £ we should use the action Tklg^i^] at k = 1/i because with this functional a tree level 
analysis is sufficient to describe the essential physics at this scale, including the relevant 
quantum effects. Hence, when we observe the space-time with a microscope of resolution 
£, we will see an average radius of curvature given by rc{£) = rdk = !/£)■ Once £ is 
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smaller than the Planck length ipi = nip^ we are in the fixed point regime where Afe oc k'^ 
so that rc{k) oc 1/k, or 

rc(£) oc e . (6.8) 

Thus, when we look at the structure of space-time with a microscope of resolution £ <^ £pi, 
the average radius of curvature which we measure is proportional to the resolution itself. 
If wc want to probe finer details and decrease £ we automatically decrease Tc and hence 
mcrcase the average curvature. Space-time seems to be more strongly curved at small 
distances than at larger ones. The scale-free relation (6.8) suggests that at distances 
below the Planck length the QEG space-time is a special kind of fractal with a self- 
similar structure. It has no intrinsic scale because in the fractal regime, i.e., when the RG 
trajectory is still close to the NGFP, the parameters which usually set the scales of the 
gravitational interaction, G and A, are not yet "frozen out". This happens only later on, 
somewhere half way between the non-Gaussian and the Gaussian fixed point, at a scale of 
the order of mp\. Below this scale, Gk and Afe stop running and, as a result, rdk) becomes 
independent of k so that rc{£) = const for £:$> £p\. In this regime {gnv)f, is fc-independent, 
indicating that the macroscopic space-time is describable by a single smooth, classical 
Riemannian manifold. 

6.2 Dynamical dimensional reduction of the graviton propagator 

As another consequence of the NGFP, ref. [13] observed that the 4-dimensional graviton 
propagator undergoes a dynamical dimensional reduction at high energies. The corre- 
sponding argument is based upon the anomalous dimension t/at = —dtlnZ^k, eq. (4.17), 
which for the RG trajectories of the Einstein-Hilbert truncation (within its domain of 
validity) have ?7iv « for ^ 0^ and r]N ~ —2 for k oo. The smooth change by two 
units occurs near k ^ mpi. 

*In the case of Type Ilia trajectories [14,97] the macroscopic fc-value is still far above fcterm, i.e., in the 
"GR regime" described in [97]. 
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This information can be used to determine the momentum dependence of the dressed 
graviton propagator for momenta S> mpj. Expanding 

Vk[g] = (167rGfe)-^ Jd^x^{-R + 2Xk} + gauge fixing (6.9) 

about fiat space and omitting the standard tensor structures we find the inverse propagator 
Okip)~^ ZN{k)p^- The conventional dressed propagator Q{p) contained in F = V^^q 
obtains from the exact Qk in the Hmit — t- 0. For > ^ rripj the actual cutoff scale is 
the physical momentum p^ itself^ so that the fc-evolution of Qk {p) stops at the threshold 
k = \f^. Therefore 

Q{p)-^ oc Zjv (/c = Vp^) p^ oc (p2)i-i (6.10) 

because ZN{k) oc k~'^ when 77 is (approximately) constant. In d dimensions, and for 
rj ^2 — d, the Fourier transform of Q{p) oc l/{p^)^~'^^^ yields the following propagator in 
position space: 

G{x;y)^- \r2T^. (6.11) 

\x-y\ 

This form of the propagator is well known from the theory of critical phenomena, for 
instance. (In the latter case it applies to large distances.) Eq. (6.11) is not valid directly 
at the NGFP. For d = 4 and 7] = —2 the dressed propagator is G{p) = l/p"^ which has the 
following representation in position space: 

Q{x;y) = -^ln{t^\x-y\) . (6.12) 

Here /x is an arbitrary constant with the dimension of a mass. Obviously (6.12) has the 
same form as a 1/p^-propagator in 2 dimensions. 

Slightly away from the NGFP, before other physical scales intervene, the propagator 
is of the familiar type (6.11) which shows that the quantity rjN has the standard interpre- 
^See section 1 of ref. [98] for a detailed discussion of "decoupling" phenomena of this kind. 
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tation of an anomalous dimension in the sense that fluctuation effects modify the decay 
properties of Q so as to correspond to a space-time of effective dimensionality 4 + rjN- 
Thus the properties of the RG trajectories imply a remarkable dimensional reduction: 

Space-time, probed by a "graviton" with <C rripj is 4-dimensional, but it appears to be 
2-dimensional for a graviton with 2> rripj. 

It is interesting to note that in d classical dimensions, where the macroscopic space- 
time is d-dimensional, the anomalous dimension at the fixed point is 77 = 2 — c?. Therefore, 
for any d, the dimensionality of the fractal as implied by t/at is d + r] = 2 [13, 15]. 

7 Fractal dimensions within QEG 

Having encountered the first evidence for the fractal nature of the effective QEG space- 
times, we now quantify their properties by computing their spectral, walk, and Hausdorff 
dimensions introduced in appendix A in the framework of QEG [49,60]. These generalized 
dimensions can readily be obtained by studying diffusion processes of test particles on the 
effective space-times. 

7.1 Diffusion processes on QEG space-times 

Since in QEG one integrates over all metrics, the central idea is to replace the classical 
return probability of a diffusing test particle Pg{T), eq. (A. 3), by its expectation value 

P(T) = (P-y(r)) = J P7PCPC'P^(r)e-^bare[7,C,C] _ (^J 

Here 7^,^ denotes the microscopic metric and iSbare is the bare action related to the UV 
fixed point, with the gauge-fixing and the pieces containing the ghosts C and C included. 
For the untraced heat kernel (A. 2), we define likewise 

K{x, x'; T) = {K^{x, x'; T)) . (7.2) 
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Following the discussion in subsection 6.1, these expectation values are most conve- 
niently calculated from the effective average action. Since defines an effective field 
theory at the scale k we know that {O^^^u)) ^ 0{{g^u)k) provided the observable O in- 
volves only momentum scales of the order of k. We apply this rule to the RHS of the 
diffusion equation, O = —A^K^{x,x';T). The subtle issue here is the correct identifica- 
tion of k. If the diffusion process involves (approximately) only a small interval of scales 
near k over which does not change much, the corresponding heat equation contains 
the operator A{k) for this specific, fixed value of k: dTK{x,x';T) = —A{k)K{x,x';T). 
Denoting the eigenvalues of A(A;o) by 6n and the corresponding eigenfunctions by this 
equation is solved by 

K{x, x'- T) = ^ 4>n{x)4>n{A exp ( - F{k^)enT) . (7.3) 

n 

Here we introduced the convenient notation F{k'^) = Afc/A^,). Knowing the propagation 
kernel, we can time-evolve any initial probability distribution p{x; 0) according to 

pix;T) = J d'^xW9o{x')K{x,x';T)p{x';0) (7.4) 

with go the determinant of {g^u)ko- If the initial distribution has an eigenfunction expan- 
sion of the form p{x; 0) = J2n Cn't>n{x) we obtain 

p{x- T) = Y, CnMx) exp ( - F{k'')£nT) . (7.5) 

n 

If the Cn's are significantly different from zero only for a single eigenvalue we are 
dealing with a single-scale problem and would identify fe^ = £n as the relevant scale at 
which the running couplings are to be evaluated. In general the (7„'s are different from 
zero over a wide range of eigenvalues. In this case we face a multiscale problem where 
different modes (j)n probe the space-time on different length scales. If A(A;o) corresponds 
to flat space, say, the eigenfunctions = are plane waves with momentum p'^, and 
they resolve structures on a length scale I of order Hence, in terms of the eigenvalue 
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6n = £p= the resolution is £ 1/ y/S^. This suggests that when the manifold is probed 
by a mode with eigenvalue £n it "sees" the metric {gij,v)k for the scale k = \fS^. Actually, 
the identification k = is correct also for curved space since, in the construction of 
Fj;., the parameter k is introduced precisely as a cutoff in the spectrum of the covariant 
Laplacian. 

As a consequence, under the spectral sum of (7.5), we must use the scale = 6n 
which depends explicitly on the resolving power of the corresponding mode. Likewise, in 
eq. (7.3), F(fc^) is to be interpreted as F{£n): 

K{x, x'; T) = J2 Mx)Mx') exp ( - F{Sn)SnT^ 

(7.6) 

(-F(A(fco))A(A;o)r)0 

n 

As in [49], we choose ko as a macroscopic scale in the classical regime, and we assume that 
at ko the cosmological constant is small, so that {g^v)ko can be approximated by the flat 
metric on W^. The eigenfunctions of A(A;o) are plane waves then and eq. (7.6) becomes 

K{x, x'; T) = J ^ e-P"^^"^^ (7.7) 

where the scalar products are performed with respect to the flat metric, {gfj,u)ko — ^tiv- 
The kernel (7.7) satisfles K(x, x';0) = 8'^{x — x') and, provided that \mip^Qp^F{p^) = 0, 
also / d'^xK{x,x';T) = 1. 

Taking the normalized trace of (7.7) within this "flat space-approximation" yields [49] 

P(r) = / ^ e-P'^^'^^ . (7.8) 



(27r) 

Introducing z = p'^, the final result for the average return probability reads 

^(^'=(S)^rW2)X ''"^'^-«p(-.FWt), (7.9) 
where F{z) = X{k'^ = z)/\kQ- In the classical case, F{z) = 1, this relation reproduces 
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the familiar result P{T) = l/(47rr)'^/^. Thus, in this case, the spectral dimension (A. 5) is 
given by Vg (T) = d independently of T. 



7.2 The spectral dimension in QEG 

We shall now discuss the spectral dimension (A. 5) for several other illustrative and im- 
portant examples. 

(A) To start with, let us evaluate the average return probability (7.9) for a simplified RG 
trajectory where the scale dependence of the cosmological constant is given by a power 
law, with the same exponent 6 for all values of k: 

X^o^k^ =^ F{z) (X //^ . (7.10) 

By rescaling the integration variable in (7.9) we see that in this case 

const 



Hence (A. 5) yields the important result 



(7.12) 



It happens to be T-independent, so that for T ^ trivially 



<i, = ^. (7.13) 



(B) Next, let us be slightly more general and assume that the power law (7.10) is valid 
only for squared momenta in a certain interval, G \z\ , Z2] , but Afe remains unspecified 
otherwise. In this case we can obtain only partial information about P{T), namely for 
T in the interval The reason is that for T € [^2""'^, z^^^] the integral in (7.9) is 

dominated by momenta for which approximately 1/p^ « T, i.e., z G [21,22]- This leads 
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us again to the formula (7.12), which now, however, is vahd only for a restricted range 
of diffusion times T; in particular the spectral dimension of interest may not be given by 
extrapolating (7.12) to T ^ 0. 

(C) Let us consider an arbitrary asymptotically safe RG trajectory so that its behavior 
for A; —> DO is controlled by the NGFP. For k > mpi the scale-dependence of the cosmo- 
logical constant is governed by the fixed point, i.e., Afe = A*A;^, independently of d. This 
corresponds to a power law with (5 = 2, which entails in the NGFP regime, i.e., for 

V,{T) = I (^NGFP regime) . (7.14) 

This dimension, again, is locally T-independent. It coincides with the T — > limit: 

ds = '^. (7.15) 

This is the result first derived in ref. [49]. As it was explained there, it is actually an exact 
consequence of Asymptotic Safety which relies solely on the existence of the NGFP and 
does not depend on the Einstein-Hilbert truncation. 

(D) Returning to the Einstein-Hilbert truncation, let us consider the piece of the Type 
Ilia RG trajectory depicted in Fig. 9 which lies inside the linear regime of the GFP. 
Newton's constant is approximately fc-independent there and the cosmological constant 
evolves according to 

Afc = Ao + J^Gofc'^. (7.16) 

Here v = {A-Kf-'^l'^id - is a scheme-dependent constant [11,14]. When k is 

not too small, so that Aq can be neglected relative to uGok'^, we are in what we shall 
call the "A;*^ regime"; it is characterized by a pure power law A^ ^ k^ with S = d. The 
physics behind this scale dependence is simple and well-known: It represents exactly the 
vacuum energy density obtained by summing up the zero point energies of all field modes 
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integrated out. For T in the range of scales pertaining to the k'^ regime we find 

Id 

'^s{T) = j^ (fc' regime). (7.17) 

Note that for every d> 2 the spectral dimension in the k'^ regime is even smaller than in 
the NGFP regime. 



7.3 The walk dimension in QEG 

In order to determine the walk dimension for the diffusion on the effective QEG space- 
times we return to eq. (7.7) for the untraced heat kernel. We restrict ourselves to a regime 
with a power law running of Afc, whence F{p^) = {LpY with some constant length-scale 
L. 

Introducing = p^r^/(^+'^) and = (x^ — x'^) /T^/^'^'^^'> we can rewrite (7.7) in the 
form 

^(-'-';^) = ^w^^(j^) (7-18) 

with the function 

d'^q 



For (5 = 0, this obviously reproduces (A.6). From the argument of $ in (7.18) we infer 
that r = \x — x'\ scales as T^l^'^+^) so that the walk dimension can be read off as 



V^{T) = 2 + 5 



(7.20) 



In analogy with the spectral dimension, we use the notation T)u,{T) rather than dw to 
indicate that it might refer to an approximate scaling law which is valid for a finite range 
of scales only. 
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For S = 0,2, and d we find in particular, for any topological dimension d. 



2 



classical regime 



4 NGFP regime 
2 + d k'^ regime 



(7.21) 



Regimes with all three walk dimensions of (7.21) can be realized along a single RG tra- 
jectory. Again, the result for the NGFP regime, = 4, is exact in the sense that it does 
not rely on the Einstein-Hilbert truncation. 

7.4 The Hausdorff dimension in QEG 

The smooth manifold underlying QEG has per se no fractal properties whatsoever. In 
particular, the volume of a d-ball B'^ covering a patch of the smooth manifold of QEG 
space-time scales as 



Thus, by comparing to eq. (A. 7), we read off that the Hausdorff dimension is strictly equal 
to the topological one: 



We emphasize that the effective QEG space-times should not be visualized as a kind 
of sponge. Their fractal-like properties have no simple geometric interpretation; they are 
not due to a "removing" of space-time points. Rather they are of an entirely dynamical 
nature, reflecting certain properties of the quantum states the system "space-time metric" 
can be in. 




(7.22) 



dn = d 



(7.23) 
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7.5 The Alexander- Orbach relation 



For standard fractals the quantities d^, d^, and dn are not independent but are related 
by [99] 

By combining eqs. (7.12), (7.20), and (7.23) we see that the same relation holds true for 
the effective QEG space-times, at least within the Einstein-Hilbert approximation and 
when the underlying RG trajectory is in a regime with power-law scaling of Afe. For every 
value of the exponent 5 we have 

The results dn = d, T)w = 2 + 5 imply that, as soon as (5 > d — 2, we have > dn and 
the random walk is recurrent then [100] . Classically {5 = 0) this condition is met only in 
low dimensions d < 2, but in the case of the QEG space-times it is always satisfied in the 
k'^ regime {6 = d), for example. So also from this perspective the QEG space-times, due 
to the specific quantum gravitational dynamics to which they owe their existence, appear 
to have a dimensionality smaller than their topological one. 

It is particularly intriguing that, in the NGFP regime, P^, = 4 independently of d. 
Hence the walk is recurrent {T>vj > dn) ior d < 4, non-recurrent for d > 4, and the 
marginal case Vyj = dn is realized if and only if d = 4, making d = 4 & distinguished 
value. Notably, there is another feature of the QEG space-times which singles out d = 4: 
it is the only dimensionality for which ^^(NGFP regime) = d/2 coincides with the effective 
dimension d^g = d + rj^ = 2 derived from the graviton propagator in section 6.2. 



8 The scale-dependence of Vg and Vw 

We now proceed by discussing the scale-dependence of the spectral and walk dimension, 
arising within the Einstein-Hilbert truncation. For this purpose, we consider an arbitrary 
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RG trajectory k i-^ {gk, Xk)- Along such an RG trajectory there might be isolated intervals 
of A;- values where the cosmological constant evolves according to a power law, Xk oc k^, for 
some constant exponents S which are not necessarily the same on different such intervals. If 

the intervals are sufficiently long, it is meaningful to ascribe a spectral and walk dimension 
to them since 5 = const implies /c-independent values Vg = 2d/ (2 + 6) and = 2 + 6. 

In between the intervals of approximately constant Vg and , where the A;-dependence 
of Afe is not a power law, the notion of a spectral or walk dimension might not be mean- 
ingful. The concept of a scale-dependent dimension Vg or V^j is to some extent arbitrary 

with respect to the way it connects the "plateaus" on which 6 = const for some extended 
period of RG time. While RG methods allow the computation of the T>s and "D^, values 
on the various plateaus, it is a matter of convention how to combine them into continuous 
functions k Vs{k),T>w{k) which interpolate between the respective values. 

8.1 The exponent 5 as a function on theory space 

In this subsection, we describe a special proposal for a /c-dependent T>s{k) and T>w{k) 
which is motivated by technical simplicity and the general insights it allows. We retain 
eqs. (7.12) and (7.20), but promote 6 — )• 5{k) to a /c-dependent quantity 

5{k) = kdkHXk). (8.1) 

When Xk satisfies a power law, oc k^ this relation reduces to the case of constant 5. 
If not, 5 has its own scale dependence, but no direct physical interpretation should be 
attributed to it. The particular definition (8.1) has the special property that it actually 
can be evaluated without first solving for the RG trajectory. The function 6{k) can be 
seen as arising from a certain scalar function on theory space, 5 = 5{g,X), whose k- 
dependence results from inserting an RG trajectory: 5{k) = S{gk,Xh)- In fact, (8.1) 
implies S{k) = kdk ln(A;^Afe) = 2 + X'^^kdkXk so that 5{k) = 2 + X'^^/3x{gk, Xk) upon using 
the RG equation kdkXk = l3\{gk,Xk)- Thus when we consider 5 as a function on theory 
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Figure 9: The g-A— theory space with the line of turning points, B, and a typical trajectory of 
Type Ilia. The arrows point in the direction of decreasing k. The big black dot indicates the 
NGFP while the smaller dots represent points at which the RG trajectory switches from increasing 
to decreasing A or vice versa. The point T is the lowest turning point, and C is a typical point 
within the classical regime. For A > 0.35, the RG flow leaves the classical regime and is no longer 
reliably captured by the Einstein-Hilbert truncation. (From [60].) 

space, coordinatized by g and A, it reads 

6{g,X) = 2 + j(3x{9,X). (8.2) 

Substituting this relation into (7.12) and (7.20), the spectral and the walk dimensions 
become functions on the ^-A-plane 

Ps(g,A) = 4^;^-?^^(g ^ P»(5,A) =4 + A-Va(<7,A). (8.3) 

To evaluate these expressions further, we use the beta functions (4.16). As we already 
discussed, the scaling regime of a NGFP has the exponent 6 = 2. From eq. (8.2) we learn 
that this value is realized at all points {g, A) where /3x = 0. The second condition for the 
NGFP, I3g = 0, is not required here, so that we have 6 = 2 along the entire line in theory 
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space: 



B = {{g,X)\Px{9A)=0}. (8.4) 



For d = 4 the curve B is shown as the bold blue line in fig. 9. Both the GFP {g, A) = (0, 0) 
and the NGFP, {g, A) = (5*, A*), are located on this curve. Furthermore, the turning points 
T of all Type Ilia trajectories are also situated on B, and the same holds for all the higher 
order turning points which occur when the trajectory spirals around the NGFP. This 
observation leads us to an important conclusion: The values 6 = 2 ■^=^ Vg = d/2, = 4 
which (without involving any truncation) arc found in the NGFP regime, actually also 
apply to all points {g, A) G B, provided the Einstein-Hilbert truncation is reliable and no 
matter is included. 



8.2 The spectral and walk dimensions along a RG trajectory 

We proceed by investigating how the spectral and walk dimension of the effective QEG 
space-times change along a given RG trajectory. Our interest is in scaling regimes where 
Vs and remain (approximately) constant for a long interval of fc-values. For the 
remainder of this subsection, we will restrict ourselves to the case c? = 4 for concreteness. 

We start by numerically solving the coupled differential equations (4.15) with the beta 
functions (4.16) evaluated with the optimized cutoff for a series of initial conditions keeping 
Ainit = ^{ko) = 0.2 fixed and successively lowering ^init = g{ko). The result is a family 
of RG trajectories where the classical regime becomes more and more pronounced. These 
solutions are substituted into (8.3), which give 'Ds{t; ginit, ^init) and 'Dyj{t; g-mit, ^init) in 
dependence of the RG-time t = In(fe) and the RG trajectory. One can verify explicitly, 
that substituting the RG trajectory into the return probability (7.9) and computing the 
spectral dimension from (A. 4) by carrying out the resulting integrals numerically gives 
rise to the same picture. 

Fig. 10 then shows the resulting spectral dimension, the walk dimension, and the lo- 
calization of the plateau-regimes on the RG trajectory in the top-left, top-right and lower 
diagram, respectively. In the top diagrams, girnt decreases by one order of magnitude for 
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Figure 10: The t = ln(A:)-dependent spectral dimension (upper left) and walk dimension (upper 
right) along illustrative solutions of the RG equations (4.15) in d = 4. The trajectories develop 
three plateaus: the classical plateau with Vs = 4,1?,^, — 2, the semi-classical plateau where Vs = 
4/3, = 6 and the NGFP plateau with Vs — 2,T>w — 4. These plateau values are indicated by 
the gray horizontal lines and connected by crossover parts. The lower figure shows the location 
of these plateaus on the RG trajectory: the classical, /c^, and NGFP regime appear between the 
points Pi and P2, P3 and P4, and above P5, respectively. (From [60].) 

each shown trajectory, starting with the highest value to the very left. As a central result, 
fig. 10 establishes that the RG flow gives rise to three plateaus where T)s{t) and T>m{t) are 
approximately constant: 

(i) For small values k, below t ~ 1.8, say, one finds a classical plateau where Ds = 4, = 2 
for a long range of A;- values. Here 5 = 0, indicating that the cosmological constant is in- 
deed constant. 

(ii) Following the RG flow towards the UV (larger values of t) one next encounters the 
semi- classical plateau where Vg = 4/3,2?^ = 6. In this case 6{k) = 4 so that Xk cx on 
the corresponding part of the RG trajectory. 

(iii) Finally, the NGFP plateau is characterized hy T>s = 2, "Dw = 4, which results from 
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the scale-dependence of the cosmological constant at the NGFP Xk oc S = 2. 

At this stage, it is worthwhile to see which parts of a typical RG trajectory realize 
the scaling relations underlying the plateau- values of Vg and Vw This is depicted in the 
lower diagram of fig. 10 where we singled out the third solution with ginit = 10~^ for 
illustrative purposes. In this case the classical plateau is bounded by the points Pi and P2 
and appears well below the turning point T, while the semi-classical plateau is situated 
between the points P3 and P4 well above the turning point. The NGFP plateau is realized 
by the piece of the RG trajectory between P5 and the NGFP. The turning point T is not 
situated in any scaling region but appears along the crossover from the classical to the 
semi-classical regime of the QEG space-times. For t < 0, the spectral dimension (walk 
dimension) increases (decreases) rapidly. In this region, the underlying RG trajectory is 
evaluated outside the classical regime at points A > 0.35. In this region of the theory 
space, the Einstein-Hilbert truncation is no longer trustworthy, so that this rapid increase 
of T>s is most likely an artifact, arising from the use of an insufficient truncation. 

The plateaus observed above become more and more extended the closer the trajec- 
tories turning point T gets to the GFP, i.e., the smaller the IR value of the cosmological 
constant. The first RG trajectory with the largest value g-mit = 0.1 does not even develop 
a classical and semi-classical plateau, so that a certain level of fine-tuning of the initial 
conditions is required in order to make these structures visible. Interestingly enough, 
when one matches the observed data against the RG trajectories of the Einstein-Hilbert 
truncation [51,97] one finds that the "RG trajectory realized by Nature" displays a very 
extreme fine-tuning of this sort. The coordinates of the turning point are approximately 
5(T ~ At ~ 10~^° and it is passed at the scale kx ~ lO"^"/??.?! « 10~^eV « (10~^mm)~^, 
so that there will be very pronounced plateau structures in this case. 
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Figure 11: Spectral dimension I?J'^'^(r) determined from random walks on a CDT space-time 
built from N = 200k simplices [56]. 

8.3 Matching the spectral dimensions of QEG and CDT 

The key advantage of the spectral dimension T>s{T) is that it may be defined and computed 
within various a priori unrelated approaches to quantum gravity. In particular, it is easily 
accessible in Monte Carlo simulations of the Causal Dynamical Triangulations (CDT) 
approach in d = 4 [54] and d = 3 [56] as well as in Euclidean Dynamical Triangulations 
(EDT) [59]. This feature allows a direct comparison between Pp^'^(T) and !D^^^(T) 
obtained within the discrete approaches and dP^'^(T) capturing the fractal properties of 
the QEG effective space-times. We conclude this section by reviewing this comparison for 
(i = 3, following ref. [60]. In particular we shall determine the specific RG trajectory of 
QEG which, we believe, underlies the numerical data obtained in [56]. 

Let us start by looking into the typical features of the spectral dimension Pp^"'"(T) 
obtained from the simulations. A prototypical data set showing T>^^'^{T) as function of 
the length of the random walk T is given in fig. 11. The resulting curve is conveniently 
split into three regimes: 

(i) For T < 20, corresponding to the left gray region in fig. 11, V'^^'^iT) undergoes rapid 
oscillations. These originate from the discrete structure of the triangulation to which the 
short random walks are particular sensitive. 

(ii) For long random walks with T > 500, the data shows an exponential fall-off. This 
^°For a similar fit in Hofava gravity see ref. [101]. 




69 



feature is due to the compact nature of the triangulation, which imphes that for long 
random walks Vf^'^{T) is governed by the lowest eigenvalue of the Laplacian on the 
compact space. This regime is marked by the right gray region in fig. 11. 
(iii) Between these two regimes, Dp^^(T) is affected neither by the discreteness nor 
the compactness of the triangulation. Since for 2?^^'^(T), determined by the flat-space 
approximation (7.8), we do not expect any of these effects to appear, we use this middle 
region to compare the T-dependent spectral dimensions arising from the two, a priori 
different, approaches. 

The comparison between the CDT-data and the QEG effective space-times is carried 
out as follows: 

(i) First, wc numerically construct a RG trajectory gk{go, Aq), Xk{go, Aq) depending on the 
initial conditions go,^o, by solving the flow equations (4.15). 

(ii) Subsequently, we evaluate the resulting spectral dimension 'D'^^'^{T; go, Aq) of the cor- 
responding effective QEG space-time. This is done by first finding the return probability 
P{T;go,Xo), eq. (7.9), for the RG trajectory under consideration and then substituting 
the resulting expression into (A. 5). Besides on the length of the random walk, the spec- 
tral dimension constructed in this way also depends on the initial conditions of the RG 
trajectory. 

(iii) Finally, we determine the RG trajectory underlying the CDT-simulations by fitting 
the parameters 50, Aq to the Monte Carlo data. The corresponding best-fit values are 
obtained via an ordinary least-square fit, minimizing the squared Euclidean distance 

500 2 
{AVsf ^ Yl {^?^''iT; 9l\ Afi*) - Vf^^iT)) , (8.5) 

T=20 

between the (continuous) function D?^*^(T; (701 Aq) and the points I?^^^(T). We thereby 
restrict ourselves to the random walks with discrete, integer length 20 < T < 500, which 
constitute the white part of flg. 11 and correspond to the regime (iii) discussed above. 

The resulting best-flt values 5'o*,Aq* for the triangulations with N = 70.000, N = 
100.000, and N = 200.000 simplices are collected in table 4. Notably, the sum over 
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70k 



0.7 X 10"^ 



7.5 X 10"^ 



0.680 



100k 



8.8 X 10^^ 



39.5 X 10-5 



0.318 



200k 



13 X 10-5 



61 X 10-5 



0.257 



Table 4: Initial conditions Jq*, Aq' for the RG trajectory providing the best fit to the Monte Carlo 
data [56]. The fit-quality {AVg)^, given by the sum of the squared residues, improves systematically 
when increasing the number of simplices in the triangulation. 

the squared residuals in the third column of the table improves systematically with an 
increasing number of simplices. By integrating the flow equation for g{k),X{k) for the 
best-fit initial conditions one furthermore observes that the points gQ^,X^^ are actually 

located on different RG trajectories. Increasing the size of the simulation N leads to a 
mild, but systematic increase of the distance between the turning point T and the GFP 
of the corresponding best-fit trajectories. 

Fig. 12 then shows the direct comparison between the spectral dimensions obtained 
by the simulations (blue curves) and the best-fit QEG trajectories (green curves) for 70k, 
100k and 200k simplices in the upper left, upper right and lower left panel, respectively. 
This data is complemented by the relative error 



for the three fits in the lower right panel. The 70k data still shows a systematic deviation 
from the classical value T)s{T) = 3 for long random walks, which is not present in the 
QEG results. This mismatch decreases systematically for larger triangulations where the 
classical regime becomes more and more pronounced. Nevertheless and most remarkably 
we find that for the 200k-triangulation that e < 1%, throughout. All three sets of residues 
thereby show a systematic oscillatory structure. These originate from tiny oscillations 
in the CDT data which are not reproduced by P?^*^(T). Such oscillations commonly 
appear in systems with discrete symmetries [70] and are thus likely to be absent in the 



vr''{T;gl\Xf) 



(8.6) 
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Figure 12: Comparison between the 3-diniensional CDT data-sets 70k (upper left), 100k (up- 
per right), and 200k (lower left) obtained in [56] (blue curves) and the best fit values for 
I?^^^(T; (7o*7 -^0*) (green curves). The relative errors for the fits to the CDT-datasets with 
N = 70.000 (circles), N = 100.000 (squares) and N = 200.000 (triangles) simplices are shown 
in the lower right. The residuals grow for very small and very large durations T of the random 
walk, consistent with discreteness effects at small distances and the compactness of the simulation 
for large values of T, respectively. The quality of the fit improves systematically for triangulations 
containing more simplices. For the N = 200k data the relative error is « 1%. (From [60].) 



continuum computation. As a curiosity, we observe that the QEG result matching the 
most extensive simulation with = 200k "overshoots" the classical value T)s{T) = 3, 
yielding Vs (T) > 3 for T > 450. At this stage, the RG trajectory is evaluated outside 
the classical regime in a region of theory space where the Einstein- Hilbert approximation 
starts to become unreliable. It is tempting to speculate that larger triangulations may 
also be sensitive to quantum gravity effects at distances beyond the classical regime. 

We conclude this section by extending V^^'~^ {T; g^^ , Aq') obtained from the 200k data 
to the region of very short random walks T < 20. The result is depicted in fig. 13 which 
displays the CDT-data V^^^{T) (blue curve), the best-fit QEG result V'^^^{T; X^^) 
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Figure 13: Extrapolation of the CDT-data for the spectral dimension measured from Monte 
Carlo simulations with 200.000 simplices [56] (blue curve) to infinitesimal random walks. The 
green and the red curve show the QEG-prediction V'^^^iT; gg', Ag') and the best-fit function (8.7) 
employed in [56], respectively. Notably, the scaling regime corresponding to the NGFP is reached 
for log(T) < —40, which is well below the distance scales probed by the Monte Carlo simulation. 

(green curve) and the best-fit function 

Vf^^-^\T) = a + be-''^, a = 3.05 , 6 = -1.02, c = 0.017, (8.7) 

employed in [56] (red curve) as a function of log(T). Similarly to the 4-dimensional case 
discussed in fig. 10, the function 2??^*^(T; (7q*, Ag*) obtained for d = 3 develops three 
plateaus where the spectral dimension is approximately constant over a long T-interval. 
For successively decreasing duration of the random walks, these plateaus correspond to 
the classical regime PP^*~^(r) = 3, the semi-classical regime where pP^*~'(T) ~ 1 and the 
NGFP regime where P?^^(T) = 3/2. The figure illustrates that P^^'^(T) probes the 
classical regime and part of the first crossover towards the semi-classical regime only. This 
is in perfect agreement with the assertion [56] that the present simulations do not yet probe 
structures below the Planck scale. The different predictions for the spectral dimensions 
obtained for infinitesimal random walks solely arise from different extrapolation schemes. 
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9 Concluding remarks 



In this article we reviewed the basic concepts underlying the gravitational Asymptotic 
Safety program and summarized evidence supporting the existence of the non-Gaussian 
renormalization group fixed point (NGFP) at the heart of this scenario. The NGFP 
constitutes a highly non-trivial feature of the gravitational renormalization group flow and 
defines a consistent and predictive quantum theory of gravity, Quantum Einstein Gravity 
or QEG for short. Based on the continuum average action approach [11], we then outlined 
our present understanding of the multifractal features characteristic for the effective QEG 
space-times. These structures are at least to some extend captured by the spectral, walk, 
and Hausdorff dimension seen by a fictitious diffusion process set up on the effective 
space-times. Their explicit computation allows a direct comparison between the results 
obtained from the continuum functional renormalization group and other approaches to 
quantum gravity like the discrete Causal Dynamical Triangulations [52-54, 57, 58] and 
Euclidean Dynamical Triangulation [59] programs. We close our review with the following 
comments: 

(1) We stress that the construction of an effective average action for gravity reviewed in 
section 3 represents a "background independent" approach to quantum gravity. Somewhat 
paradoxically, this background independence is achieved by means of the background 
field formalism: One fixes an arbitrary background, quantizes the fiuctuation field in this 
background, and afterwards adjusts g/^i, in such a way that the expectation value of the 
fluctuation vanishes: hfj,i, = 0. In this way the background gets fixed dynamically. The 
results of refs. [28,29] suggest that a proper implementation of "background independence" 
might be crucial for Asymptotic Safety. 

(2) The combination of the effective average action with the background field method has 
been successfully tested within conventional field theory. In QED and Yang-Mills type 
gauge theories it reproduces the known results and extends them into the non-perturbative 
domain [39,41]. 

(3) The coexistence of Asymptotic Safety and perturbative non-renormalizability is well 
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understood. In particular upon fixing gf^^, = rj^v and expanding the trace on its RHS in 
powers of G the FRGE reproduces the divergences of perturbation theory; see ref. [5] for 
a detailed discussion of this point. 

(4) It is to be emphasized that in the average action framework the RG flow, i.e., the 
vector field is completely determined once a theory space is fixed. As a consequence, 
the choice of theory space determines the set of fixed points V* at which asymptotically 
safe theories can be defined. Therefore, in the Asymptotic Safety scenario the bare action 
S related to V* is a prediction of the theory rather than an ad hoc postulate as usually 
in quantum field theory. Ambiguities could arise only if there is more than one suitable 
NGFR 

(5) According to the results available to date, the Einstein-Hilbert action of classical 
General Relativity seems not to play any distinguished role in the Asymptotic Safety 
context, at least not at the conceptual level. The only known NGFP on the theory space 
of QEG has the structure V* = Einstein-Hilbert action + "more" where "more" stands for 
both local and non-local corrections. So it seems that the Einstein-Hilbert action is only 
an approximation to the true fixed point action, albeit an approximation which was found 
to be rather reliable for many purposes. 

(6) A solution to the FRGE, as such, does not define a regularized functional integral; a 
priori it is not even clear whether or not a path integral representation exists for a given 
RG trajectory. In [77] it was shown how, upon specification of a UV regularized measure, 
the information contained in r^_^(-x3 may be used to obtain the cutofl^ dependence of the 
bare action. A thorough understanding of this relationship is essential in comparing the 
average action to other approaches, such as Monte Carlo simulations of discrete models of 
gravity, for instance. 

(7) Any quantum theory of gravity must reproduce the successes of classical General 
Relativity. As for QEG, it cannot be expected that this will happen for all RG trajectories 
in S\jY, but it should happen for some or at least one of them. Within the Einstein-Hilbert 
truncation it has been shown [97] that there actually do exist trajectories (of Type Ilia) 
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which have an extended classical regime and are consistent with all observations. 

(8) In the classical regime mentioned above the dynamics of space-time geometry is un- 
affected by quantum effects to a very good approximation. In this regime the familiar 
methods of quantum field theory in curved classical space-times apply, and it is clear 
therefore that effects such as Hawking radiation or cosmological particle production are 
reproduced by the general framework of QEG with matter. 

(9) Coupling free massless matter fields to gravity, it turned out [21] that the fixed point 
continues to exist under very weak conditions concerning the number of various types of 
matter fields (scalars, fermions, etc.). No fine tuning with respect to the matter multi- 
plets is necessary. In particular Asymptotic Safety does not seem to require any special 
constraints or symmetries among the matter fields such as supersymmetry, for instance. 

(10) Since the NGFP seems to exist already in pure gravity it is likely that a widespread 
prejudice about gravity may be incorrect: its quantization seems not to require any kind 
of unification with the other fundamental interactions. 

(11) To the best of our knowledge, the effective space-times arising within QEG carry 
a multifractal structure [60]. Some of these features are captured by the spectral, walk, 
and Hausdorff dimension seen by a fictitious diffusion process set up on the effective 
space-time. The resulting Hausdorff dimension is constant and equal to the topological 
dimension of the (background) space-time. Thus the fractal properties do not originate 
from the QEG space-times "loosing points" at short distances but rather represent a 
genuine dynamical effect of quantum field theory. In contrast to the Hausdorff dimension 
the spectral dimension and the walk dimension seen by the diffusion process depend on 
the diffusion time. Fig. 10 identifies three regimes in which these generalized dimensions 
are constant for a wide range of scales, the classical, semi-classical and NGFP regime, 
which are connected by short crossovers. 

(12) In section 8.3 we performed a direct comparison between the spectral dimension of 
the 3-dimensional effective QEG space-times with the one measured in Causal Dynamical 
Triangulations (CDT) [56]. Notably, the best-fit RG trajectory reproduces the CDT data 
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with approximately 1% accuracy for the range of diffusion times where the simulation data 
is reliable. Notably current Monte Carlo data neither probes the semi-classical plateau 
nor the scaling regime of the NGFP. Determining the spectral dimension for infinitesimal 

random walks by extrapolating the leading quantum corrections to the classical regime 
may therefore miss the imprints of the NGFP on the structure of space-time. 
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A Generalized dimensions on classical manifolds 

Investigating random walks and diffusion processes on fractals, one is led to introduce 
various notions of fractal dimensions, such as the spectral or walk dimension [100]. These 
notions also prove useful when characterizing properties of space-time in quantum gravity, 
and we will review these concepts in the remainder of this section. 

A.l The spectral dimension 

To start with, consider the diffusion process where a spin-less test particle performs a 
Brownian random walk on an ordinary Riemannian manifold with a fixed classical metric 
g^j,i,{x). It is described by the heat kernel Kg{x,x';T) which gives the probability density 
for a transition of the particle from x to x' during the fictitious time T. It satisfies the 
heat equation 

dTKg{x,x';T) = -AgKg{x,x'-T), (A.l) 
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where Ag = —D"^ denotes the Laplace operator. In flat space, this equation is easily solved 

by 

K,{x,x'-T) = J ^e^^'-(— ')e-P^^ (A.2) 

In general, the heat kernel is a matrix element of the operator exp(— TAg). In the random 
walk picture its trace per unit volume, 

PgiT) = j d'^x^f^Kg{x,x-T) = IV exp(-TAg) , (A.3) 

has the interpretation of an average return probability. Here V = j d'^X\/g{x) denotes the 
total volume. It is well known that Pg possesses an asymptotic early time expansion (for 
T ^ 0) of the form Pg{T) = (AttT)-'^/^ Y.n=o ^nT"", with An denoting the Seeley-DeWitt 
coefficients. From this expansion one can motivate the definition of the spectral dimension 
ds as the T-independent logarithmic derivative 



d, = -2 



dlnPg{T) 



dlnT 



(A.4) 

T=0 



On smooth manifolds, where the early time expansion of Pg{T) is valid, the spectral 
dimension agrees with the topological dimension d of the manifold. 

Given Pg{T), it is natural to define a, in general T-dependent, generalization of the 
spectral dimension by 

According to (A.4), we recover the true spectral dimension of the space-time by considering 
the shortest possible random walks, i.e., by taking the limit ds = limT^o T^s{T). Note that 
in view of a possible comparison with other (discrete) approaches to quantum gravity the 
generalized, scale-dependent version (A.5) plays a central role. 
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A. 2 The walk dimension 



Regular Brownian motion in flat space has the celebrated property that the random 
walker's average square displacement increases linearly with time: (r^) oc T. Indeed, 
performing the integral (A.2) we obtain the familiar probability density 

K{x,x'-T) = (47rr)-'^/2exp (-%?^) (A.6) 

with a{x,x') = ^\x — x'S^ half the squared geodesic distance between the points x^x' . 
Using (A.6) yields the expectation value (r^) = (.x^) = J d^xx^ K{x^{^\T) oc T. 

Many diffusion processes of physical interest (such as diffusion on fractals) are anoma- 
lous in the sense that this linear relationship is generalized to a power law (r^) oc T^/'^"' 
with dy, ^ 2. The interpretation of the so-called walk dimension d^ is as follows. The 
trail left by the random walker is a random object, which is interesting in its own right. 
It has the properties of a fractal, even in the "classical" case when the walk takes place 
on a regular manifold. The quantity is precisely the fractal dimension of this trail. 
Diffusion processes are called regular li dy, = 2, and anomalous when dy, ^2. 



A. 3 The Hausdorff dimension 

Finally, we introduce the Hausdorff dimension dH- Instead of working with its mathe- 
matically rigorous definition in terms of the Hausdorff measure and all possible covers of 
the metric space under consideration, the present, simplified definition may suffice for our 
present purposes. On a smooth set, the scaling law for the volume V{r) of a d-dimensional 
ball of radius r takes the form 

V{r) oc r'^" . (A.7) 
The Hausdorff dimension is then obtained in the limit of infinitely small radius, 

dH ^ hm . (A.8) 

r-i^o Inr 
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Contrary to the spectral or walk dimension whose definitions are linked to dynamical 
diffusion processes on space-time, there is no such dynamics associated with dn- 
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